AD-A238  012 

■■111111 

RL-TR-91  -92 
Final  Technical  Report 
July  1991 


DESIGN  AND  ANALYSIS  OF  MULTI¬ 
SENSOR  SEQUENTIAL  DETECTION 
SYSTEM 


University  of  Virginia 

Quan-Ming  Chen  and  Demetrios  Kazakos 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


91-04686 


% ELECTE 
^JUL  15 


Rome  Laboratory 
Air  Force  Systems  Command 
Griffiss  Air  Force  Base,  NY  1 3441  -5700 

91 


1  A 

l. 


This  report  has  been  reviewed  by  the  Rome  Laboratory  Public  Affairs  Office 
(PA)  and  is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RL-TR-91-92  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


VINCENT  C.  VANNICOLA 
Project  Engineer 


APPROVED: 


JAMES  W.  YOUNGBERG,  LtCol,  USAF 
Deputy  Director  of  Surveillance 


FOR  THE  COMMANDER: 

IGOR  G.  PLONISCH 

Directorate  of  Plans  &  Programs 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  Rome  Laboratory 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  please 
notify  RL(OCTS  )  Griffiss  AFB  NY  13441-5700.  This  will  assist  us  in  maintaining  a 
current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


orm  Approved 
OMB  No.  0704-0188 


Putafc:  'acortng  bLfOTi  far  tn»  coMcncn  of  rforrfon  a  «grrno  to  w« g»  i  hOLf  par  r capons*  noucfriQ  ty»  tm»  for  r»»»ng  nanjcoon*  saarcnng  amstng  aaa  sojces 
gafwng  and  nm  taring  an*  oata  n— am  ana  canMMig  ana  raviawng  tfucotacocn  d  rformam  Send  cammarts  rtgadng  tna  Oli dan  ascmaa  or  any  octaer  aspect  or  :ms 
colacDon  c*  rtormam  ndudng  -  ^XT***-”* for  raducing  tn»  tuOarv  to  Waarr^on  Haackquartara  Sarvcaa,  Oraacroa  far  rformaoon  Oparncna  ana  Reports.  1 21 5  „erterson 
Davw  Higfsrray.  SuU  1 204.  Aifr^nn,  VA  22202-4302.  and  to  thi  Oflta  of  Managaman  and  0uOg*.  PapanaoK  Rarljctpn  Proiag  (0704-01 88),  Wag-ngorv  PC  20503 


1.  AGENCY  USE  ONLY  (Leave  Blank) 


2  REPORT  DATE 

July  1991 


4.  TITLE  AND  SUBTITLE 

DESIGN  AND  ANALYSIS  OF  MULTI-SENSOR  SEQUENTIAL 
DETECTION  SYSTEM 
6.  AUTHOR  (S)  ~ 

Quan-Ming  Chen,  Demetrios  Kazakos 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
University  of  Virginia 
Department  of  Electrical  Engineering 
Charlottesville  VA  22901 


9.  SPONSORING/MONTT ORING  AGENCY  NAME(S)  AND  ADDRESSES) 

Rome  Laboratory  (0CTS) 

Griffiss  AFB  NY  13441-5700 


a  REPORT  TYPE  AND  DATES  COVERED 
Final  Jan  89  -  Dec  89 


a  FUNDING  NUMBERS 

C  -  F30602-81-C-0169 
PE  -  62702F 

-  PR  -  4506 

TA  -  00 
WU  -  TR 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


RL-TR-91-92 


1 1 .  SUPPLEMENTARY  NOTES 

Rome  Laboratory  Project  Engineer:  Vincent  C.  Vannicola/OCTS/ (315)  330-4437 


1 2a  DISTRIBUTION/AVAILABILfTY  STATEMENT  1 2b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


1  a  ABSTRACT  (Moran  200  wort*) 

A  multiple  sensor  system  is  considered  under  binary  hypothesis  environments.  All 
sensors  are  assumed  independent,  and  the  observed  data  is  also  independent.  Re¬ 
ceived  data  is  quantized  and  then  sent  to  the  fusion  center  to  determine  whether 
a  target  is  present.  The  Sequential  Probability  Ratio  Test  is  employed  in  the 
fusion  center.  The  objective  is  to  find  an  optimal  system  by  minimizing  the  expected 
number  of  observations.  Both  two-level  and  four-level  quantizer  are  used  in  the 
process  of  finding  the  optimal  system.  Numerical  evaluations  are  made  to  find  the 
quantizer  which  minimize  the  expected  number  of  observations  that  are  required  to 
decide  the  presence  of  the  target.  System  simulations  are  also  performed  to  confirm 
the  results. 


1 4.  SUBJECT  TERMS 

Detection,  Fusion,  Communications,  Radar 


15  NUMBER  OF  PAGES 

120 _ 

14  PRICE  CODE 


1 7.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


NSN  7540-01 .280-5500 


1  a  SECURITY  CLASSIFICATION  1 9.  SECURITY  CLASSIFICATION  20.  UMITATION  OF  ABSTRACT 

°F  LJNCLA^llFIED  °WfIED  U/L 


Stvoara  Form  298  rRev  2  891 
Prmat* d  by  ANSI  Std  Z 39  1 8 
298*1 02 


Abstract 


A  multiple  sensor  system  is  considered  under  binary  hypothesis  environments. 
All  sensors  are  assumed  independent,  and  the  observed  data  is  also  independent. 
Received  data  is  quantized  and  then  send  to  the  fusion  center  to  determine  whether  a 
target  is  present.  The  Sequential  Probability  Ratio  Test  is  employed  in  the  fusion 
center.  The  objective  to  find  an  optimal  system  by  minimizing  the  expected  number 
of  observations.  Both  two-level  and  four-level  quantizer  are  used  in  the  process  of 
finding  the  optimal  system.  Numerical  evaluations  are  made  to  find  the  quantizer 
which  minimize  the  expected  number  of  observations  that  are  required  to  decide  the 
presence  of  the  target.  System  simulations  are  also  performed  to  confirm  the  results. 


i 


Accession  For 

NT IS  GFA&I 
OTIC  TAB 

Unannounced 

Justif icatioi 

□ 

□ 

"j- - 

Dlstr 

Aval 

Dlst 

r 

ifcut ion 

labillt 

Avail 

Spec 

/ 

f  Codes 
md/d*'" 

Lai 

Table  of  Contmts 


Chapter  1.  Introduction .  1 

1.1  Literature  Review  and  Goals .  1 

1.2  Overview  of  Chapters .  2 

1.3  Environment .  3 

1.4  Assumption .  3 

1.5  Advantages  of  the  Sequential  Test .  4 

Chapter  2.  General  Analysis  of  a  Two-Sensor  System .  5 

2.1  The  Model  and  Configuration  .  5 

2.1.1  Decision-Making  at  the  Individual  Sensor .  5 

2.1.2  Decision-Making  at  the  Fusion  Center  .  5 

2.1.3  .  8 

2.2  .  10 

2.2.1  General  Definition  of  SPRT . •. .  10 

2.2.2  Calculation  of  the  Thresholds  .  10 

2.3  SPRT  for  Two  Distributed  Sensors  .  1 1 

2.3.1  Definition  of  the  LR .  1 1 

2.3.2  SPRT .  12 

2.3.3  Calculation  of  the  Thresholds  A'  and  B' .  12 

i  i 


2.3.4  Derive  the  Expected  Value  of  N  ( Ej  {  N  }) .  15 

2.4  Evaluate  E;  {N}  Under  Quantization .  16 

2.4.1  Definition  of  the  G  Function .  18 

Chapter  3.  Analysis  of  a  Two-Sensor  System  with  Two-level  Quantizer .  19 

3.1  The  Model  and  Configuration .  19 

3.1.1  The  Model  of  Sensor .  19 

3.2  Evaluate  the  Expected  Value  of  N  (Ej  {N}) .  21 

3.2.1  An  Alternative  to  Derive  Ej  {N } .  22 

3.2.2  The  G  Function .  24 

3 3  Minimization  of  G  Under  the  Double  Exponential  Assumption .  24 

3.3.1  Probability  Density  Function .  25 

3.3.2  The  Sequential  Probability  Ratio  Test .  25 

33.3  Calculation  on  G  function .  26 

3.4  Numerical  Evaluation  of  G .  29 

Chapter  4.  Analysis  of  a  Two-Sensor  System  with  Four-level  Quantizer  .  37 

4.1  The  Model  and  Configuration .  37 

4.1.1  The  Model  of  Sensor .  37 

4.2  Evaluate  the  Expected  Value  of  N  (Ej  {N}) .  39 

4.2.1  The  G  Function .  40 

4.3  Minimization  of  G  Under  the  Double  Exponential  Assumption .  41 

iii 


4.3.1  Probability  Density  Function .  41 

43.2  The  Sequential  Probability  Ratio  Test .  41 

43.3  Calculation  on  G  function .  42 

4.4  Numerical  Evaluation  of  the  G  Function .  54 

Chapter  5.  Comparison  of  the  Performance  of  Sequential  Systems .  62 

5.1  Performance  of  Individual  System .  62 

5.1.1  G  of  Two-level  Quantizer  System .  62 

5.1.2  G  of  the  Four-level  Quantizer  System .  64 

5 2  Comparison  of  System .  65 

Chapter  6.  System  Simulations .  68 

6.1  Simulation  Method .  68 

68 

6.1.1  Derivation  of  The  Double  Exponential  Environment . 

6.1.2  Simulation  and  Discussion .  69 

Chapter  7.  Conclusion  .  76 

Appendix  A.  Program  for  Evaluation  of  G  Function  of  Two-level  Quantizer  System  .  78 

Appendix  B.  Program  for  Evaluation  of  G  Function  of  Four-level  Quantizer  System . 

Appendix  C.  Program  for  Evaluation  of  G  Function  of  Two-level  Quantizer  System  .  98 

References  105 

i  v 


List  of  Figures 


Figure  2.1  Model  of  m-level  Quantizer .  6 

Figure  2.2  The  m-level  Quantized  Probability  Space  of  Sensor  1  .  7 

Figure  2.4  Model  of  the  Sequential  Detection  System  .  9 

Figure  3.1  Quantized  Probability  Space  of  Two-level  Quantizer  System 

.  20 

Figure  3.2  G  with  a  =  P  =  0.01  .  31 

Figure  3.3  G  with  a  =  p  =  0.05 .  32 

Figure  3.4  G  with  a  =  .01 ,  p  =  .05 .  3  3 

Figure  3.5  G  with  a  =  .01,(3  =  .05  and  means  at  -2,1  and  -3,2  .  34 

Figure  3.6  G  with  a  =  .05, (3  =  .01  and  means  at  -3,3  and  -5,5  .  35 

Figure  3.7  G  with  a=  .01, P  =  .05  and  means  at  -3,3  and  -5,5  .  36 

Figure  4.1  Four-level  Quantized  Probability  Space .  38 

Figure  4.2  G  of  FQSS  with  a=p=. 01  .  56 

Figure  4.3  G  of  FQSS  with  Various  T i  and  T2 . : .  5 1 

Figure  4.4  G  of  FQSS  with  a=P=.05  .  58 

Figure  4.5  G  of  FQSS  with  a=.01  P=.05  .  59 

Figure  4.6  G  of  FQSS  with  =  .01  =  .05  and  60 

.  means  at  -2.1  and  3.2 

Figure  4.7  G  of  FQSS  with  Ti=-.1T2=-.2  and  means  at  -3,3  and  -5,5 .  61 

Figure  6.1  N  of  Two-level  Quantizer  System .  70 


V 


Figure  6.2  N  of  Four-level  Quantizer  System .  *  71 

Figure  63  Enlarged  N  with  a  =  J3  =  0.01 .  72 

Figure  6.4  Enlarged  N  with  a  =  0.01  (J  =  0.05 .  73 

Figure  63  Enlarged  N  with  a  =  0.05  (S  =  0.05 .  74 

Figure  6.6  Enlarged  N  with  a  =  .05  p  =  .01  and  means  at  -3 3  and  -S3 .  75 


vi 


List  of  Tables 


Table  5.1  Tabulated  Data  of  Minimum  G  with  Two-level  Quantizer  . 

Table  5.2  Minimum  G  with  Two-level  Quantizer  a  =  0.05  (5  =  0.01  .  64 

Table  5.3  Tabulated  Data  of  Minimum  G  with  Four-level  Quantizer .  65 

Table  5.4  Minimum  G  with  Four-level  Quantizer,  a  =  0.05  P  =  0.01  .  65 

Table  5.5G  with  Various  Ti  and  T2 .  66 

Table  5.6  Improvement  on  G  .  67 


vii 


CHAPTER  1 


Introduction 


1.1.  Literature  Review  and  Goals 

The  Sequential  Probability  Ratio  Test  (SPRT)  was  first  developed  by  Wald.  In 
[1],  he  presented  the  general  theory  of  the  sequential  analysis  and  SPRT.  Over  the 
years,  there  are  many  books  have  been  written  about  the  sequential  analysis.  In  par¬ 
ticular,  Dr.  Siegmund  has  discussed  this  topic  in  a  very  precise  and  condensed  form 
in  his  book  [2]. 

In  [3],  sequential  detection  based  on  simple  quantization  has  been  analyzed.  In 
this  paper,  the  classical  ruin  problem  in  probability  theory  was  applied  to  analyze  the 
sequential  dead-zone  limiter  detector  and  the  sequential  four-level  detector.  How¬ 
ever,  the  author  did  not  study  the  the  effects  of  the  quantization  level  of  a  quantizer 
on  the  sequential  system,  and  his  results  were  only  applied  to  a  single  sensor.  Chair, 
Hoballah  and  Varshney  [4]  applied  the  sequential  detection  theory  to  decentralized 
system  where  multisensors  were  used.  A  global  decision  was  made  based  on  the 
local  decision  of  each  sensor.  Local  decision  rules  were  given  by  the  likelihood  ratio 
test  (LRT),  and  the  authors  used  the  Neyman-Pearson  approach  to  derive  the  optimal 
rules  for  each  detector. 

In  our  research,  we  also  considered  the  problem  of  the  distributed  sequential 
detection  system.  However,  instead  of  using  LRT  at  each  sensor  as  mentioned  in 
Chair,  Hoballah,  and  Varshney’s  paper,  we  simply  quantize  the  observations  into  m- 
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level  and  sent  them  to  the  fusion  center,  the  fusion  center  employs  a  sequential  pro¬ 
cess  that  has  the  option  to  make  a  final  decision  on  whether  or  not  to  continue  the 
process  by  taking  one  more  observation  from  each  sensor.  Our  goals  were  to  obtain 
the  analytical  expression  on  the  expected  number  of  samples  that  are  required  to 
make  a  decision  in  terms  of  the  probabilities  associated  with  the  quantized  observa¬ 
tions,  and  to  investigate  the  effects  of  different  quantizers  on  the  expected  number. 
We  are  particularly  interested  in  the  systems  which  have  two-level  and  four-level 
quantizers.  The  performances  are  compared  through  numerical  evaluations  and  sys¬ 
tem  simulations. 

12.  Overview  of  Chapters 

Chapter  1  contains  a  literature  review  in  which  some  of  the  important  papers  are 
discussed.  It  also  describes  the  goals  of  the  study,  and  the  general  overview  of  each 
chapter  is  given.  The  conditions  and  general  assumptions  are  also  mentioned. 

A  discussion  on  a  two-sensor-system  with  m-level  quantizer  is  presented  in 
chapter  2.  The  sequential  probability  ratio  test  is  stated  for  the  system.  The  expres¬ 
sion  of  the  expected  number  of  observations  is  also  determined,  and  G  function  is 
defined  in  terms  of  expected  number  of  observations. 

In  chapter  3,  the  emphasis  is  on  the  two-sensor  system  with  two-level  quantizer. 
By  assuming  that  the  distribution  function  of  the  observations  has  a  double  exponen¬ 
tial  character,  the  numerical  evaluations  are  performed  on  the  expected  number  of  the 
system  for  various  quantizers,  and  the  optimal  systems  are  found,  the  effects  of  the 
parameters  of  the  double  exponential  on  the  system  are  also  considered.  The  results 
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are  plotted,  and  the  evaluation  program  is  in  appendix  A. 

Work  similar  to  that  done  in  the  previous  chapter  is  done  for  the  system  with 
four-level  quantizer  in  chapter  4.  The  analytical  expression  as  well  as  the  numerical 
evaluation  of  the  expected  number  are  performed  for  different  four-level  quantizers 
and  parameters  of  the  double  exponential  distribution  function.  The  data  is  plotted  in 
the  end  of  the  chapter,  and  the  program  of  this  study  is  in  appendix  B. 

The  results  from  chapter  3  and  chapter  4  are  presented  in  chapter  5.  Each 
system’s  results  are  compared  to  the  other  results.  The  detailed  discussion  of  each 
system  is  also  given  in  this  chapter. 

In  chapter  6,  the  optimal  quantizers  that  are  found  in  previous  chapters  are  used 
to  simulate  the  two-sensor  system.  The  results  are  listed  and  plotted.  In  this  chapter, 
the  derivation  of  the  random  environment  of  double  exponential  is  also  discussed. 
The  results  arc  shown  in  the  end  of  the  chapter,  and  the  simulation  program  is  in 
appendix  C. 

In  chapter  7,  final  comments  and  conclusions  are  made  on  the  study. 

1.3.  Environment 

The  environment  consists  two  hypothesis  Hi  and  Ho-  shows  that  a  target  is 
present,  and  Hq  indicates  that  no  target  is  present. 

1.4.  Assumption 

We  assume  that  the  data  observed  at  a  sensor  is  independent  and  identical  distri¬ 
buted.  It  is  also  assumed  that  the  observations  of  one  sensor  are  independent  from 


the  others. 


1.5.  Advantages  of  the  Sequential  Test 

The  sequential  test  is  considered  optimal  in  the  sense  that  it  minimizes  the 
expected  sample  size  both  under  Hi  and  Ho  among  all  tests  have  no  larger  error  pro¬ 
bability  with  independent  and  identical  distributed  observation  [2].  Its  average  sam¬ 
ple  size  is  smaller  than  the  fixed  sample  size  required  by  the  Neyman-Pearson  test  for 
the  same  performance. 
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CHAPTER  2 


General  Analysis  of  a  Two-Sensor  System 


2.1.  The  Model  and  Configuration 

2.1.1.  Decision-Making  at  the  Individual  Sensor 

There  arc  two  hypotheses  Hx  and  H<j.  At  each  sensor,  the  hypotheses  are 
presented  by  the  probability  density  functions.  The  sensor  quantizes  the  incoming 
data  into  m-levels.  The  m-level  quantizer  shown  in  Figure  2.1  is  precalculated 
according  to  the  optimal  rule  which  will  be  derived  in  later  chapters.  As  shown  in 
Figure  2.2,  each  quantized  level  is  associated  with  two  conditional  probabilities.  If 
Xj  falls  between  a*  and  ak_lt  the  probabilities  associate  with  X,  are  pKJ  and  pK0- 
Figure  2.3  shows  the  quantized  regions  under  hypotheses  Hj  and  Hq  of  sensor  two. 
The  total  number  of  bits  that  is  transmitted  to  the  fusion  center  depends  on  the 
number  of  the  quantized  level.  For  a  four-level  quantizer,  it  requires  two  bits  of  data. 

2.1.2.  Decision-Making  at  the  Fusion  Center 

At  the  fusion  center.  Sequential  Probability  Ratio  Test  (SPRT)  is  employed. 
The  decision  is  made  by  comparing  the  two  predetermined  thresholdes  A  and  B 
where  A  is  less  than  B.  The  fusion  center  will  request  more  information  from  the 
sensors  if  it  cannot  decide  whether  the  target  is  present. 
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The  m-level  Quantizer 


Figure  2.1  Model  of  m-level  Quantizer 


Quantized  Regions  for  Sensor  1 


Figure  2.2  The  m-level  Quantized  Probability  Space  of  Sensor  1 


Quantized  Region  for  Sensor  2 
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Figure  2.3  The  m-level  Quantized  Probability  Space  of  Sensor  2 


2.1  J.  The  Overall  Process  of  the  Model 

As  mentioned  in  chapter  1,  sensors  are  independent  each  other.  The  incoming 
signals  are  first  quantized  into  m-levels,  and  then  send  to  the  fusion  center.  Depend¬ 
ing  on  whether  the  likelihood  ratio(LR)  at  the  fusion  center  is  greater  than  B  or  less 
than  A,  the  final  decision  is  made.  A  target  is  either  detected  if  LR  is  greater  than  B, 
or  the  target  is  not  present  if  LR  is  less  than  A.  If  the  LR  is  laid  in  between  the  thres¬ 
holds  A  and  B,  the  center  requests  each  sensor  to  send  one  more  quantized  observa¬ 
tion  as  shown  in  Figure  2.4. 
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Figure  2.4  Model  of  the  Sequential  Detection  System 


2.2.  Definition  of  the  Sequential  Probability  Ratio  Test 


2.2.1.  General  Definition  of  SPRT 

If  xi,X2,...  is  a  sequence  of  random  variables  with  join  density  function  fn,  and 
there  are  two  hypotheses  Hj  and  Hq  such  that: 

Ho  •  fn  ~  fbn(xl»x2»  •  •  *  »xn) 

Hj  .  fn  fln(xl*x2»  •  •  •  »xn) 

ln  is  defined  as  likelihood  ratio: 

ln  =ln(xl*x2»  •  •  •  >xn) 

_  fln(xl  »x2»  •  •  •  »xn) 

fbn(xl  »x2»  •  •  •  » xn) 

Let  the  thresholds  A  and  B  are  chosen  such  that  0  <A  <  B  <  <»,  and  start  observing 
data  XJ.X2, ....  sequentially  until  the  random  time  N  which  is  the  first  n  such  that 
either  I„  >  B,  or  1„  <  A.  If  ln  £  B,  then  Hi  is  selected,  and  if  ln  £  A,  then  Ho  is 
selected.  This  can  be  expressed  mathematically  as  following: 

first  n  £  1  such  that  ln  £  (A,B) 

oo  if  ln  e  (A,  B)  for  all  n  >  1  <2-2-4) 

Stop  sampling  at  the  time  N  and  if  N  <  «> 

Reject  H!  if  ln  >  B 
Accept  Hq  if  ln  <  A 

2.2.2.  Calculation  of  the  Thresholds 

Thresholds  A  and  B  were  calculated  in  terms  of  the  error  probabilities  of  first 
and  second  kind  (a  and  P,  respectively)  in  [1], 


(2.2.1) 

(2.2.2) 


(2.2.3) 
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For  given  a  and  {3  where 

a  =  PoOn  -  B) 
P  =  Pj{lN  <  A) 

A  and  B  can  be  expressed  as  following: 


a 


(2.2.5) 

(2.2.6) 

(2.2.7) 

(2.2.8) 


23.  SPRT  for  Two  Distributed  Sensors 

23.1.  Definition  of  the  LR 

Extending  the  definition  of  the  LR  in  (2.2.3)  to  two  sensors,  the  LR  of  this  sys¬ 
tem  can  be  written  as: 

,  fln(xl»x2 . xn;yi.Y2 . Yn)  . 

In  =  7-7 - - - : - 7  (2.3.1) 

M)n(xl*x2»  •  •  •  »xn*  yi»Y2»  •  •  •  »yn) 

where  fin(xltx2,  •  •  •  »xn'»  Yi.Y2.  •  •  •  *Ya)  and  f<)n(xi.x2 . xn;  Yi.Y2>  •  •  •  »Yn)  arc 

join  density  function  of  observed  random  sequence  xi,x2, . . .  ,xn  from  sensor  1  and 
yi,y2, . . .  ,yn  from  sensor  2  under  hypotheses  Hj  and  Hq  respectively  (see  Figure 
2.3). 

Applyng  the  assumption  that  the  observations  of  one  sensor  are  independent 
from  the  other,  the  above  equation  (2.3.1)  can  be  re-written  as: 

.  fln(xl»x2>  •  •  •  > xn)^ln(Yl *Y2 . Yn) 

ln  =  7—7 - rr— - r  (2.3.2) 

*On(xl  *x2»  •  •  •  *  xn)^On(y  1  *y2*  •  •  •  *  y^) 

Furthermore,  if  the  individual  observation  of  each  sensor  is  independently  distributed, 
(2.3.2)  becomes: 
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_  "  flnWflnfrk) 

"  jLl  f0n(Xk)f0n(yk) 

Taking  log  on  both  sides  to  arrive  at  the  log  likelihood  ratio: 


n 


login  =  I  log 

k=l 


fln(*k)fln(yk) 

f0n(xk)f0n(yk) 


232.  SPRT 


(2.3.3) 


(2.3.4) 


Using  the  LR  that  is  defined  in  the  previous  section,  the  SPRT  can  be  stated  as 
follows: 


Sampling  the  random  sequence  xi,X2,...  and  yi,y2»—  sequentially  until  the  ran¬ 
dom  time  N,  such  that  login  is  greater  or  equal  to  b,  or  less  than  and  equal  to  a.  The 
hypothesis  is  accepted  if  login  >  b  and  the  hypothesis  Hq  is  accepted  if  login  S  a, 
where 


a  =  logA' 
b  =  logB' 

A'  and  B'  will  be  derived  in  the  next  section. 


233.  Calculation  of  the  Thresholds  A'  and  B' 

As  defined  in  equation  (2.2.5),  a  is  the  error  probability  of  first  kind,  and 
a  =  Prob{  decide  Hj  I  Hq  ) 

om 

=  £Prob{  stop  at  n  and  decide  H!  I  Ho) 

n=l 


=  £P0(N  =  n,  in>B0 

n=l 

(2.3.6) 

Let  Bn  be  the  subset 

of 

n-dimensional 

space 

where 

A'  <  lk(*i .  •  •  • ,  Xkiyi , . . . ,  yk)  <  B' 

for 

every 

k  <  n. 

and 
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Hence, 


the 


set 


ln(xl »  •  •  •  >  Xn,y  1 »  •  •  •  >  Yn)  —  B  • 

{  N  =  n,  ln  >  B' }  =  {(xlt . . .  ,xn;ylf . 


••.yn)6  Bn}. 


Let  An  be 

the  subset 

of 

n-dimensional 

space 

where 

A  < ln(xj , . . . ,  xk,y i , . . 

.  .,yk)<B' 

for 

every 

k  <  n 

and 

ln(x1,...,xn;yi,...,yn)<A/.  Hence  the  set 

{  N  =  n,  ln  <  A  }  —  { (xj , . . . ,  xn,yj , . . . , yn)  s  An ) . 

Now  the  definition  in  equation  (2.3.6) 


a=  £  J  fon(xi,...,xn,yi . yn)dx!  •  •  •  dxndyi  •  •  •  dyn 

n=l  X}«  •  •  •  •  •  •  •  •  .yne®n 

oo 

=  Z  J fOn(xi ,  •  •  • ,  x„)fo„(yi ,  •  •  • ,  yn)dxi  •  •  •  dx„dyi  •  •  •  dyn.  (2.3.7) 

n=l  Bn 

Using  the  condition  that  1„  >  B', 


1  °° 

a  <  —  Z  j  fin(xi»  •  •  •  .Xn)f|„(yi,  •  •  •  ,yn)dxi '  *  *  d*ndyi  •  •  •  dyn 
1  00 

=  — Z  lfin(xi . xn,y1,...,yn)dxi  •••dxndy1  -  --dyn 

B  n=lBn 

=  -Sr£Pl(N=".ln2B') 

B  n=l 

=  -^-Prob{  decide  Hj  I  Hj ) 

B 


=  —  (1  -  Prob{  decide  Hq  I  Hi  }) 

B' 

1-B 

=  (2.3.8) 


By  the  same  method, 
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P  =  Prob{  decide  Ho  I  Hi } 

=  £Prob{  stop  at  n  and  decide  Ho  I  Hj } 

n=l 

=  £p1(N  =  n.ln<A') 

n=l 

oo 

=  x  J  fin(xi,...,xn,yi,...,yn)dxi  •••dXndyi  • 

n=l  xi . Xn.yj . yneAn 

0O 

=  X  J  fin(xi, .  •  •  ,xn)fin(yi, .  •  •  ,yn)dxi  *  •  •  dXndyi  •  •  •  dyn. 

n=lAn 

oo 

^  A'X  S fon(xi, . . .  ,xn)f0n(yi, . . . , yn)dxi  •  •  •  dx„dyj  •  •  •  dyn 

n=lA„ 

=  A'£  J fon(xi, . . .  ,xn,ylf . . .  ,yn)dx1  •  •  •  dx„dyi  •  •  •  dyn 

n=lAn 

=  A'£p0{N  =  n,ln<A') 

n=l 

=  A'Prob{  decide  Hq  I  Hq} 


"dyn 


=  A'(l  -  Prob{  decide  Hj  I  Ho)) 

=  (1  -  a)A\  (2.3.9) 

The  a  and  p  can  be  approximated.  Since  the  equations  (2.3.8)  and  (2.3.9)  have  ine¬ 
qualities  only  because  ln  does  not  have  to  reach  the  boundaries  A'  and  B'  exactly 
when  it  first  leave  (A',B')>  we  can  ignore  this  discrepency  and  treat  (2.3.8)  and 
(2.3.9)  as  equalities.  Then 


1-p 
a=  — —■ 
B' 


(2.3.10) 
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(3  =  A'(l  -  a) 


(2.3.11) 


2i3.4.  Derive  the  Expected  Value  of  N  (  E,  {  N  }) 

To  derive  Hi  {  N  } ,  it  is  necessary  to  look  at  the  equation  (2.3.4).  By  taking  the 
expected  value  on  both  sides  of  the  equation,  it  becomes: 


Eilog  ln  —  Ef 


Vi  fln(xk)fln(yk) 
>1  °8  f0n(xk)f0n(yk)_ 


where  Ej  {  }  is  the  conditional  expectation  under  hypothses  H;,  and  i  =  0, 1. 


(2.4.1) 


Wald’s  identity  [1]  said  that  if  Vi»V2,—  is  i-i.d  with  mean  value  |i  =  EVi,Let  M 

be  any  integer  valued  random  variable  such  that  (M  =  n)  is  an  event  determined  by 

conditions  on  Vi,V2,  ...,Vn  for  all  n=l,2,...,  and  assume  that  EM  <  then 
M 

E(£Vk)  =  pEM. 

k=l 


From  (2.4.1),  we  can  define  Vk  as: 


Vk  =  log 


fln(Xk)fln(yk) 

f0n(Xk)f0n(yk) 


(2.4.2) 


Since  fin(xk)>fin(yk)»f0n(xk).  and  fon(yk)  are  all  i.i.d,  log*1";  is  also  i.i.d. 

f0n(xk)fon(yk) 

Using  the  definition  of  Wald’s  identity  on  the  right  side  of  the  equation  (2.4.1),  it 


gives 


Ejlog  ln  =  mHi  { N }  =  EjViEi {  N  )  (2.4.3) 

On  the  other  hand,  the  left  side  of  the  equation  (2.4.1)  can  also  be  approximated  as  a 

two-valued  random  variable  taking  on  the  values  a  and  b  by  the  same  reasoning  used 

to  derive  (2.3.10)  and  (2.3.11). 
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Eilog  ln  =  aPiGn  <  A')  +  bPi  (2.4.4) 

After  applying  equations  (2.3.10)  and  (2.3.11),  and  doing  some  simplifications,  the 

equation  (2.4.4)  can  be  re-written  as: 


E0logln  =  a(l-a)  +  b^— & 

B 


=a(1_j^r)+b 


1-A' 


B'-l 
Br  —  A' 
B' 


_  a(B7  —  1)  +  b(l  -  A') 
B'  -  A' 


and 


E1Logln  =  aA/(l-a)  +  b(l-p) 

_  aA'^'-D  +  bB'd-A') 

B'-A' 

By  substitution,  equation  (2.4.3)  can  be  written  as: 

Ein=e;‘(Vii^M1 

B  -A 


E0N  =  E51{V1} 


a(B' -  1)4- b(l  -  A7) 
(B'-A') 


2.4.  Evaluate  F.j  {N}  Under  Quantization 


(2.4.5) 


(2.4.6) 


(2.4.7) 

(2.4.8) 


As  mentioned  in  section  2.1,  the  system  described  here  utilizes  a  quantizer  at 
each  sensor  to  quantize  the  received  signals.  Assuming  a  m-level  quantizer  is  used, 
thus  the  space  is  divided  into  m  intervals.  At  the  jth  interval  of  sensor  1,  there  are 
two  conditional  probabilities  associated  with  it,  namely  P(j  I  Ho)  and  P(j  I  Hi).  The 
general  conditional  probability  of  both  sensors  are  defined  below: 
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For  sensor  1: 

na 

X 

II 

X 

II 

X3 

(2.5.1) 

P(X,  =j  1  Ho)  =  pjo 

(2.5.2) 

For  sensor  2: 

P(yi  =1  1  HO  =  qn 

(2.5.3) 

P(yi=llHo)  =  qio 

(2.5.4) 

where  1,  j  =  Using  these  definitions  in  (2.5.1),  (2.5.2),  (2.5.3),  and  (2.5.4), 

the  equation  (2.4.2)  for  Vj  can  be  re-written  as  follows: 


V!  =  log 


PjiQii 

PjoQio 


(2.5.5) 


Thus  Ej  {Vj }  becomes 


Ei{V1}=Ei{log-^-} 

PjOfilO 

= ZXP(  *1  =  j.  yi  =  n  HiJiog^1- 
j  l  PjOfilO 

=  ILP(xi  =j  I  hi)P(yi  =1 1  HOlog-^- 

j  i  Pjoqio 

j  i  PjoRio 

where  i  =  1,0;  j,le(l,...,m)  by  in\oking  appropriate  independence  assumptions 


(2.5.6) 


If  the  above  equation  is  substituted  into  (2.4.7)  and  (2.4.8),  the  conditional 
expectation. of  N  can  be  expressed  as: 


,,  XT  aA'(B'-  l)  +  bB'(l  -A'),™  ,  Pjiqji, 

Ei N  =  — - — - (IZPjiqn 

B  -  A  j  i  PjoQio 


,,  xr  a(B'  -  1)  +  b(l  -  A')  ln  PjiqiK 

e0n  = - — — — - (LZPjoqioiog - ) 


(B'-A') 


j  l 


Pjoqio 


(2.5.7) 

(2.5.8) 
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2.4.1.  Definition  of  the  G  Function 


In  the  previous  section,  EjN  and  EqN  have  been  evaluated.  In  order  to  find 
optimal  solution  for  the  entire  system,  we  first  defined  a  function  G, 


G  =  £EiN. 

i=0, 1 

Now  we  express  the  G  function  in  terms  of  Ej  N  and  EqN: 


a(B'-l)  +  b(l-A'),™  n  .  Pjiqn  . 

G  + - ^ — 7^ - (XZPjoqiolog— — ) 

(B  -  A)  j  i  Pjoqio 


(2.5.9) 


If  a  and  (3  are  given,  then  A'  and  B'  are  known  from  equations  (2.3.8)  and 
(2.3.9).  The  G  function  only  depend  on  the  probability  of  quantization  levels.  Thus, 
to  have  an  optimal  system,  one  must  to  find  the  optimal  quantizers  for  both  sensors 
that  will  minimize  the  equation  G. 


In  next  chapter,  we  will  study  the  effect  of  the  two-level  quantizer  on  the  sys¬ 
tem. 
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CHAPTER  3 


Analysis  of  a  Two-Sensor  System  with  Two-level  Quantizer 


3.1.  The  Model  and  Configuration 

The  system  that  is  analyzed  here  is  closely  related  to  the  one  in  chapter  two. 
However,  in  chapter  two,  the  signals  are  quantized  into  m-level  at  each  sensor.  The 
main  concern  in  this  chapter  is  to  study  the  two-level  quantizer  system  both  theoreti¬ 
cally  and  numerically. 

3.1.1.  The  Model  of  Sensor 

Since  sensor  1  and  sensor  2  are  identically  process  the  incoming  data,  it  is  rea¬ 
sonable  to  choose  an  arbitrary  sensor  as  an  example.  Taking  a  look  at  sensor  1, 
shown  in  the  figure  3.1,  a  two-level  quantizer  with  thresholds  Ti  divides  both  proba¬ 
bility  functions  fn( )  and  foi  ( )  into  two  regions.  The  received  signal  Xk  is  quantized 
to  1  if  it  falls  below  the  threshold  Tj,  and  it  is  quantized  to  2  otherwise.  For  the  kth 
signal,  the  conditional  probabilities  associate  with  each  quantized  level  are  listed 


below: 

For  sensor  1 : 

p2i  =P(Xk  is  above  the  threshold  Tj  1  =  P(xk  =  2  I  Hj)  (3.1.1) 

p  n  =  P(  Xk  is  below  the  threshold  T j  IH!)  =  P(xk  =  l  IHj)  (3.1.2) 

P20  =  P( Xk  is  above  the  threshold  Ti  I  Ho)  =  P(xk  =  2  I  Hq)  (3. 1 .3) 

p  io  =  P(  Xk  is  below  the  threshold  T!  I  H o)  =  P(xk  =  1  I  Ho)  (3.1.4) 
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For  sensor  2: 
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(3.1.5) 


q2i  =  P(  Yk  is  above  the  threshold  T2  I  Hj )  =  P(yk  =  2  I  Hj ) 

qi !  =  P(  Yk  is  below  the  threshold  T2  1  Hj )  =  P(yk  =  1  I  Hj )  (3.1 .6) 

Q20  -  P(  Yk  is  above  the  threshold  T2  I  H0)  =  P(yk  =  2  I  Ho)  (3.1.7) 

qio  =  P(  Yk  is  below  the  threshold  T2  i  Ho)  =  P(yk  =  1  I  Ho)  (3.1.8) 

3.2.  Evaluate  the  Expected  Value  of  N  (E;  {N}) 

The  general  expression  of  E;  {N}  has  been  evaluated  in  chapter  2. 

To  derive  the  equation  of  Ei(N}  for  the  system  with  two-level  quantizer,  one  must 
find  Ej{Vi)  first.  Using  the  probabilities  defined  in  the  last  section,  the  equation 
(2.5.6)  can  be  expanded  as: 


Ei{v1}  =  2:xpjiqmog 

j  1 


Pjiqn 

pjoqio 


=  Piiquiog 


Pnqn 

Pioqio 


+ 


Pnq2ilog 


Pnq2i 

pioq2o 


where 


+ 


Paqiiiog 


P2iqn 

P20q20 


+  P2iq2ilOg 


P21  q2i 
P20q20 


Pli  =  P  li » 

P2i  =P2il 

qii  =qii'» 

q2i  =  qk 

and  i  =  0,1.  After  some  simplifications, 


Ei  { Vj }  =  P21  log 


P20P11 

P21P10 


+  q2ilog 


020^11 
<121  tJlO 


+  log 


pioqio 

Pnqn 


P20P11 

E0[V1}=p20loge-r- 

P21P10 


+  q2olog 


qzoqn 

q2i  qio 


+  log 


Pioqio 

Pnqn 


(3.2.1) 


(3.2.2) 

(3.2.3) 
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Then  the  Ej  (N)  can  be  expressed  as  follows: 


E!(N}  = 


aA'(B'  - 1)  +  b(B'-A/) 


log 


Ptoqio 


B'-A' 
l-i 


,  P20P11  ,  q2oqn 

P2ilog-  — +  q21Iog- 


P21P10 


q2iqio 


E0{N}  = 


Pnqn 

a(B'-  1)  +  b(l  -  AO 


B'-A' 
'i-i 


.  P20P11  .  q2oqn 

P2oiog— _  +  q2oiog- 


P21P10 


q2iqio 


+  log 


Pioqio 


Pnqn 


(3.2.4) 


(3.2.5) 


3.2.1.  An  Alternative  to  Derive  Ei  (N) 


In  the  case  of  two-level  quantizer,  if  the  signal  is  quantized  into  -1  or  1  depend¬ 
ing  on  whether  it  falls  below  or  above  the  threshold  Tj ,  then  the  signal  is  an  antipodal 
signal.  Let’s  define 


Sn  =  xi  +  x2  +  ' '  ’  +  xn  for  sensor  1  (3.2.6) 

Sn  =  yi  +  y2  +  •  •  •  +  yn  for  sensor  2  (3.2.7) 

Notice  here  that  {Sn}  is  a  case  of  random  walk.  If  assuming  a*  is  the  number  of  l’s 

among  n  signals,  then  n  -  a,  is  the  number  of  -l’s. 


Sn  =  ai  -(n-aO 
Sn  =  a2  -  (n  -  a2) 

Now  expressing  a;  and  n  -  a*  in  terms  of  n,  Sn,  and  Sn: 
n  +  Sn 

ai  =  — 

n-Sn 

n-aj  =  — — 


(3.2.8) 

(3.2.9) 


(3.2.10) 

(3.2.11) 


22 


(3.2.12) 


a2  = 


n  +  Sn 
2 


n  -  a2 


n-  Sn 
2 


(3.2.13) 


In  the  section  3.1,  it  was  mentioned  that  an  incoming  signal  is  quantized  into  1 
and  2.  However,  no  matter  what  value  we  assign  to  the  quantization  level,  the  proba¬ 
bility  that  associates  with  it  is  the  same.  Thus 


P(Xl=  1  1  Ho }  =  p20 

(3.2.14) 

P{Xl=-l  1  Ho }  =  Pio 

(3.2.15) 

P{Xl=  1  IH1)=p21 

(3.2.16) 

P{X!  =-l  1  Hi)  =Pll 

(3.2.17) 

P{yi  =  1  1  Ho)  =q2o 

(3.2.18) 

P{yi  =-l  i  Ho}  =  qio 

(3.2.19) 

P{yi=  1  IH1)=q21 

(3.2.20) 

P{yi=-1  IH1}=qu. 

(3.2.21) 

The  LR  becomes 


j  ft  (^i  »x2, .  •  ■ ,  xn)fj  (y i  ,y2,  •  •  • ,  Yn) 

"  fo(Xi,X2 . Xn)fo(yi,y2.---.)'n) 

[P(X|  =  1  IHi)],1[P(x,=-l  1  H|)l"— 1  [P(yi  =  1  I  H,  ll^Ply!  =-l  I  Hi)]1"” 

“  [P(x,  =  l  IHo)]*'[P(xi=-l  IHo)r*'[P{yi  =  l  IHo)]n[P{yi=-l  IHoir” 

n+Sn  n-Sn  tH-Sn  n~Sn 

[P2l]  2  [Pill  2  [Q2l3  2  [qill  2  . . 

= — ^ ^ <3-2-22> 

[P20]  2  [Pio]  2  [q2o]  2  Cqio]  2 

Taking  logrithm  on  both  sides  of  the  equation  (3.2.22), 
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,  ,  n  +  Sn ,  P21  .  n-Sn,  P11  .  n  +  Sn  q2i  n-s„  qn 

login  =  — - — log—  +  — — log—  +  — — log - +  — — log - 

2  P20  2  P10  2  q20  2  qio 

And  using  the  definition  of  expectation  and  Wald’s  identity, 


1  P21  Pi  1 

Hi  {logln }  =  —  [(Ein  +  EjnEiX^log - +  (Ein-EinEix1)log-^— 

2  P20  P10 


+  (Ein  +  EinEjyOlog—  +  (Esn  -  EinEiy^log— 
020  Qio 


(3.2.23) 

Since  E1x1=2p2i~l,  E1y1=2q2i-1,  E0x1=2p2o-1,  and  E0yi=2q20-1, 
Ei  (N)  and  Eq{N}  can  be  simplified  to  exactly  the  equations  (3.2.4)  and  (3.2.5). 


This  gives  us  an  interesting  thought  that  there  is  more  than  one  approach  to  a 
problem.  In  this  case,  we  can  think  that  the  two-level  quantization  system  can  be 
considered  as  a  random  walk  problem. 


3.2.2.  The  G  Function 

Utilizing  the  equations  derived  in  the  last  section,  the  G  function  is  expressed  by 
simple  substitution. 


G  =  E1N  +  E0N 

,  aA'(B'-l)  +  bB'(l-A')  ,  ,  a(B'-  l)  +  b(l  -  A') 

-El  tv,}  —  +Eo  (V.)  —  (3.3.1) 

where  A'  and  B'  are  defined  in  chapter  2,  and  Ej  {V} )  and  Eo{Vj }  arc  shown  in 


(3.2.4)  and  (3.2.5). 


3J.  Minimization  of  G  Under  the  Double  Exponential  Assumption 

By  applying  the  double  exponential  distribution  to  the  G  function,  one  can  show 
how  the  system  performs  numerically  and  what  the  optimal  system  will  be.  The  rea¬ 
son  that  a  double  exponential  distribution  function  is  chosen  is  that  it  is  easy  to 
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manipulate. 


3.3.1.  Probability  Density  Function 
For  sensor  1: 


fj(X)  =  ^-e  XllX  R12' 

Ho:  f0(X)  =  ^-e-XllX-^11 


(3.5.1) 

(3.5.2) 


For  sensor  2: 


fj  (Y)  =  -^-e  Xz  1  Y~^22 1 
H0:  f0(Y)  =  ^-e‘X2lY‘Wl1 


(3.5.3) 

(3.5.4) 


3.3.2.  The  Sequential  Probability  Ratio  Test 


The  log  LR  is  defined  in  (2.3.4).  By  applying  it  to  the  two-level  quantizer  case. 


n 


1°g1n=  ZIog 
k=l 


fi (xk)fi(yk) 

fo(xic)fo(yk) 


=  £iog 


k=l 


pjiqn 

Pjoqio 


where  j,l=  1,2. 


(3.5.5) 


Now,  the  two  random  sequence  xi,X2,—  and  yi,y2,—  are  sampled  sequentially 
until  the  random  time  N  where  N  is  the  first  n  that  log  ln  is  not  between  a  and  b.  At 
the  time  N, 


Accept  Hj  if  log  ln  >  b. 
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Accept  Hq  if  log  ln  <  a. 

Since  a  and  J3  are  given,  and  they  are  defined  in  (2.3.6)  and  (2.3.9),  then 


a  =  logA' 

=  1°gT^-  (3.5.6) 

b  =  logB' 

=  log^,  (3.5.7) 

a 

333.  Calculation  on  G  function 


In  order  to  calculate  the  function  G,  it  is  easier  to  evaluate  Ei  {N}  and  Eq{N} 
first.  From  (3.2.4)  and  (3.2.5),  Et  {N}  and  Ej  (N)  are  derived.  By  applying  the  dou¬ 
ble  exponential  functions. 


E,{N}  = 


aA'(B'  -  1)  +  bKB'  -  A') 


B'  -A' 


.  P20P11  ,  q20Qn  ,  PioQio 

P2ilOg -  +  q2ll0g  _  _  ■  +  log 


aA'(B'~  l)  +  b(B'  -  A') 


B'-A' 


P21P10 


^-1  — X.1 1  x  — H-12 1 


tfciQio  PiiQn 


-1 


r  1 

i^-e 


log— 

Ot 


^  ll  1 
j  A|  c~Xt i x~nii  1  j^Le-Xii*-n,2i 

T|  2  2 


|  1  X-H121  f  F-Xi  I  X-|1|1 1 

p.  2 


Jre 


J 

+  f  ~e"Xl1  y-11221  log— 

j  /  00 


-x2l  y-H2l  1  f  ^  -*-2 1  y-^22' 

2  i.  2 


t2 


J  C~X2[  y-U22 1  J  ^e-*2l  y-wi  I 


T2 


Tl  X  t2  y 

J  A’1  c~Xi  I  x-mi  I  J  ^2  C~X2 1  y - H2i 1 


+  log 


Tl  *  T2  . 

J  _Le"Xl  1  *-tM2l  C  y-^221 


I  2 


-1 


(3.2.8) 
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a(B  -  l)  +  b(l  -  A  ) 


I  -Xi  I  x-iin  I 


M  I 


a(B  -  1)  +  b(l  -  A  )  |  f  *1  -Xi  i  *-|in  i 

l  e 

1  j  ^Le~x i  1  x-wa'  J  _2±e-xi 1  *-mi 1 

Ti  ^  2 


‘^2.  ~>-2ly-H2ll 


.  f  ^2  —X2 1  y  —  M-21  *  1 

+  J  Te  lo§~ 

l 


^2-e_x2 1  y-u22'  r  ^le-X2i  y-Mi 1 

9  J  9 


+  log: 


T»  1  T2 

j  _Le-xi 1  x-wi 1  j  _e-x2 1  y-wi  i 


!  \  ‘2  * 

_^1  -Xl!x-Hl2l  f  4,2  .->-2 1  y -422 1 


(3.2.9) 


Then  G  =  Et  {N}  +  Eo(N}  can  be  calculated.  Below  are  the  list  of  calculated 
probabilities  of  all  quantized  levels  for  sensor  1: 

If  Ti  >  Hi2, 


....  f  ^Le~Xl  I  x-^12  I  _  _1  -Xl(Tj  -4i2) 


P2' =  J ' Te 


■r 


1 i  i 

_  f  ' 1  c~xl 1  x-4121  —  |  I  c~xl(Ti -412) 

Jo  9 


f  —Xt  I  x-411 1  1  -Xi(Tt  —mi) 

P20=J— e  =  -t  * 


_  r  ^Le-Xi  i  x-mi  i  _  j  _  _l_e-x1(Ti-4n) 
Jo  o c 


If  Hu  <Ti  £  Hi2» 


_  _  i  1  -M(Tl-Hl2) 

pii-i-ye 

_  _  1  -^1(Ti-H12) 

P11  -  ye 

n  _  1  -^l(Ti-nu) 

P20  -  ye 

n  _  1  _  1  -M(Tl-mi) 

P10  -  1  ye 

IfTi  <Hn, 

_  _  1  1  -*-l(Ti  -m2) 

P21  =  1  -  ye 

_  _  _L  -^l(Tl-|ll2) 

PH  -  ye 

_  _  1  _  l-MTl-Ull) 

P20  -  1  ~  ye 

_  _  -mi) 

P10  -  2 

For  sensor  2,  the  probabilities  of  quantized  level  are  similar  to  those  in  sensor  1. 
If  T2  >1122. 


_  f  c— *2  I  y-^22  1  _  J_e~*-2(T2-^22) 


_  _  r  ^2  --Vji  y-M2*  _  1  J__-^-2('i'2 —^22) 

qn  -  J  ~^~e  -  1  ~  ye 


n  _  f  -X2I  y-H21  I  _  1  -*2(T2-H2l) 
020  -  ~  Te 


_  r  ^2  c~*-2 1  y-H21  1  _  1  _  J_g-^2(T2-H2l) 

Jo  0 


If  JJ-21  <T2  —  P22> 
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_  _  t  1  -X2(T2~H22) 

q2i  -  1  -  ye 

qn  =  ye‘l2(T2‘K2) 

_  _  1  -*-2(T2-H2l) 

q20  "  T 

n  _  1  _  *  -*2(T2-«l) 

qio  - 1  ye 

If  T2  ^  fl2l » 


n  -  1  J_^2(T2-^22) 

<l2i  -  1  -  ye 

_  1  -X2(T2-H22) 
qu  "  2 

_  _  i  1  -^2(T2-ll2l) 

020  -  1  ~  ye 
n  _  1  -X2(T2-Jl2l) 

qio  -  ye 


3.4.  Numerical  Evaluation  of  G 


The  optimal  sequential  system  that  is  studied  here  is  defined  as  having  the 
minimum  amount  of  data  needed  to  determine  whether  the  target  is  present  under  the 
hypotheses  and  Hq.  Since  G  is  the  sum  of  the  conditional  expectation  of  the  total 
number  of  observation  N,  the  expected  N  under  both  hypotheses  will  be  minimized  if 
the  minimum  of  G  is  found.  Thus  the  purpose  of  the  numerical  evaluation  is  to  find 
the  best  quantizers  for  both  sensors  which  minimizes  the  G  function. 

During  the  evaluation,  different  parameters  such  as  |ii2,  Jin,  |i2i>  H22>  <*.  and  P 

are  used  to  determine  the  effect  on  the  system.  The  values  of  Xj  and  X2  remain  the 
« 

same  throughout  the  study,  and  Xi  =  0.25  and  X2  =  0.125.  By  varing  Tj  and  T2,  the 
minimum  G  is  found. 
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Figure  3.2  shows  the  G  function  with  respect  to  T2  and  Tj.  The  means  used 
here  are  symmetric  about  the  y-axis,  and  (in  =  -1,  \i\2  =  1,  fJ.21  =  -2,  and  (122  =  2. 
The  error  probabilities  a  =  p  =  0.01 

As  it  shown,  the  G  function  is  at  its  minimum  when  the  both  Tj  and  T2  of  the  quan¬ 
tizers  are  equal  to  zero. 

If  the  means  of  the  density  functions  is  kept  the  same,  and  the  error  probabilities 
increase  such  that  a  =  P  =  0.05,  the  G  is  shown  in  Figure  3.3.  The  plot  in  figure  3.4 
gives  the  relation  between  G,  Tj,  and  T2  with  a  =  0.01  and  p  =  0.05,  while  the  means 
are  still  as  same  as  in  the  previouse  figures. 

In  Figure  3.5,  the  characteristic  of  G  with  means  that  are  not  symmetrical  about 
the  y-axis  is  studied.  The  means  of  the  density  functions  are 

(in  =  -2,  P12  =  1,  (i2i  =  -3,  II22  =  2,  and  the  error  probabilities  a  =  0.01  and 
P  =  0.05.  From  this  figure,  it  is  clear  that  the  optimal  quantizer  has  its  quantized  level 
selected  such  that  it  is  in  the  center  of  the  mean. 

Figure  3.6  gives  another  perspective  when  the  (in  =  -3,  (ii2  =  3, 

(i2i  =  -5,  P22  =  5,  a  =  0.05,  and  P  =  0.01.  The  optimal  quantized  level  is  not  in  thf 
center  of  means,  T!  =-0.1,  T2  =  -0.2.  Again  in  Figure  3.7  shows  the  same  conclu¬ 
sion  with  a  different  set  of  error  probabilities,  a  =  0.01,  P  -  0.05.  The  G  is  minimum 
when  Ti  =0.3  and  T2  =0.1. 
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CHAPTER  4 


Analysis  of  a  Two-Sensor  System  with  Four-level  Quantizer 


4.1.  The  Model  and  Configuration 

The  process  of  this  sequential  detection  system  is  the  same  as  in  chapter  2  and 
chapter  3.  The  two  independent  sensors  send  quantized  observations  to  the  fusion 
center  which  employs  a  sequential  detection  scheme  to  determine  whether  a  signal  is 
present  or  not. 

The  only  difference  is  that  the  sensors  use  four-level  quantizers  instead  of  two- 
level  ones.  The  expected  number  of  data  will  be  reduced,  since  the  data  from  the 
four-level  quantizer  contain  more  information  than  the  data  from  the  two-level  quan¬ 
tizer. 

4.1.1.  The  Model  of  Sensor 

As  shown  in  Figure  4.1,  each  quantizer  consists  of  two  variables,  namely  Tj  and 
D,  where  i  =  1,2.  It  divides  the  probability  space  into  four  regions.  These  are 
(-<»,  Tj-Dj],  (Ti-Dj,T;3,  (Ti,Ti+Di],  and  (Ti  +  Dj,+°°).  After  the  sensor 
receives  a  signal,  depending  on  which  region  it  falls  in,  the  received  data  is  quantized 
to  1, 2,  3,  or  4.  Then  this  quantized  value  is  sent  to  the  fusion  center  for  processing. 

Below  are  the  list  of  probabilities  associated  with  each  quantized  level: 

For  sensor  1 : 
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Four-level  Quantizer  System 


Sensor  2 


P4i  =P(X1e(T1  +Dlf-H«)  IH1)  =  P(x1  =4  I  HO  (4.1.1) 

p31  =P(X1e(T1,T1+D1]  1  H1)  =  P(x1  =3  IH^  (4.1.2) 
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P2i  =P(X1e(T1  -Di.Tj]  I  H1)  =  P(x1  =2  I  H,) 
Pll  =P(X1€(-«>, Ti  -DO  I  H1)  =  P(x1  =  1  I  HO 
p40  =  P( X,  €  (Tj  +  Di ,  4oc)  |  Ho)  =  P(Xl  =  4  I  Ho) 
p30  =  P(X1e(T1,T1 +DO  IHq)  =  P(Xl  =  3  I  Hq) 
p20  =  P(  X! €  (T!  -  Dj ,  TO  I  Ho)  =  P(xj  =  2  I  Ho) 
pl0=P(X1e(-o,T1-D1]  I  Ho)  =  P(xj  =  1  I  H0) 

For  sensor  2: 

q4i  =  P(  Yl€ (Tj  +  D2,  +oo)  I  HO  =  P(X!  =  4  I  HO 
q31=P(Y1€(T1,T1+D2]  IH0  =  P(X!=3  IHO 
q21  =P(Y1€(T1  — D2, TO  !  H0  =  P(X!  =2  I  HO 
q„=P(Y1€(^»,T1-D2]  IH0  =  P(X!  =  1  IHO 
q40  =  P(  Yx  e  (Ti  +  D2,  -Ko)  I  Ho)  =  P(yi  =4  I  Ho) 
qso  =  P(  Yj  €  (T! ,  Tj  +  D2]  I  Ho)  =  P(yi  =  3  I  H0) 
q20  =P(  Yj6 (Ti  - D2 ,  T:]  I  Hq)  =  P(yi  =  2  I  Ho) 
q10  =  P(  Yj  e  (-~,  -  D2]  I  Ho)  =  P(yi  =  1  I  Hq) 


4.2.  Evaluate  the  Expected  Value  of  N  (Ej  {N}) 

Using  the  general  expression  of  Ej{VO  in  (2.5.6),  one 
the  four-level  quantizer  system. 


4  4 


Pjiqn 


Ei{vo  =  ziPjiqiilog-  „ 

j=n=i  Pjoqio 

.  Pnqn  ,  .  Pnq2i  , 

=  Piiqnlog— —  +  Piiq2ilog— — - —  +  Pnq3ilog 


Pioqio 


Pioq2o 


(4.1.3) 

(4.1.4) 

(4.1.5) 

(4.1.6) 

(4.1.7) 

(4.1.8) 

(4.1.9) 

(4.1.10) 

(4.1.11) 

(4.1.12) 

(4.1.13) 

(4.1.14) 

(4.1.15) 

(4.1.16) 


can  derive  Ei{VO  for 


Pnq3i 

pioq3o 
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.  PllQ41  .  P2iqil  ,  P21021 

+  Pliq4ilOg - +  P2iqiil°g -  +  P2iq2ilog — — 

Pl0q40  P20q20  P20q20 

.  P2iq3i  .  P2iq4i  ,  ,  P3iqn 

+  P2iq3ilog - +  P2iq4ilog  +  p3iqiilog—— — 

P20q30  P20q40  P30qi0 

.  P3iq2i  ,  P3iq3i  ,  P3iq4i 

+  P3iq2ilog - +  p3iq3ilog-  -  +  P3iq4iIog—  ----- 

P30q20  P30q30  P30q40 

,  P4iqn  ,  P4iq2i  ,  ,  P4iq3i 

+  p4iqiilOg - +  P4iq2ilOg - +  P4iq3iiOg - 

P40q20  P40q20  P40q30 


+  P4iq4i^Og 


p4iq4i 

P40q40 


Then  the  Ej  { N }  can  be  expressed  as  follows: 


.  {N)  .  aA-(B--l)  +  b(B--A-)[^pji^i,og^,]— 
B  -  A  |j=ii=i  Pjoqio 

_  a(B'-l)  +  b(l-A')fii  , 

Eo<N)  =  - — dr~r. - -  ISPjoqioiog-r— 

B  —  A  Lpii=i  Pjoqio 


where  a,  b,  A',  and  B'  are  all  defined  in  chapter  2. 


4.2.1.  The  G  Function 


(4.2.1) 


(4.2.2) 


(4.2.3) 


In  chapter  2,  the  G  function  is  defined  in  (2.6.1).  By  appling  the  functions 
Ei  (N)  and  Eq  {N}  derived  in  the  last  section,  G  becomes 


G  =  E1N  +  E0N 

,  aA'(B'  -  1)  +  bB'(l  -  A')  .  a(B'  -  1)  +  b(l  -  A') 

=  E‘  (V>) - i - +  E°  fVl)  B'-A' 

aA'(B'  -  1)  +  b(B'  -  A')  1*4, 4,  _  _  Pji^il”1 
= - — 77 -  IlPjiqnlog— — - 

B  -  A  1=11=!  Pjoqio 


.  a(B'  -  l)  +  b(l  -  A')  Pjiqii 

+ - 3 — 77 -  EZPjoqiolog— — 

B  -  A  j=n=i  Pjoqio 


(4.3.1) 
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4J.  Minimization  of  G  Under  the  Double  Exponential  Assumption 


4.3.1.  Probability  Density  Function 

Double  exponential  density  functions  are  applied  to  the  equation  G  to  find  the 
optimal  system.  We  assume  that  the  double  exponential  probability  density  functions 
are  as  same  as  in  chapter  3. 


For  sensor  1: 


Hit  fj (X)  =  -^-e  1  x-^12 1 

Ho: 


(4.4.1) 


(4.4.2) 


For  sensor  2: 


Hi:  fjL  (Y)  =  -^-e-*2 1  Y_W2 1 

Hq:  f0(Y)  =  ^-e-Xj,Y-«" 


(4.4.3) 


(4.4.4) 


4.3.2.  The  Sequential  Probability  Ratio  Test 


For  four-level  quantizer,  the  LR  is  defined  as: 


,  ,  fi(xk)fi(yk) 

logl”=Jilogw^7 

=  £  iog-Si^ 
k=l  Pjoqio 


(4.4.5) 


where  j,l=  1,2, 3, 4.  Then  the  SPRT  can  be  restated  as  sampling  the  two  random 
sequences  X!,x2,...  and  yj,y2,...  sequentially  until  the  random  time  N,  where 
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N  =  first  n  >  1  such  that  loglng  (a,b) 

=  «>  if  log  ln  6  (a,  b)  for  all  n. 

At  the  time  N, 


Accept  Hj  if  logln  >  b, 

Accept  Ho  if  logln  <  a. 

4 33.  Calculation  on  G  function 

Equations  (4.2.2)  and  (4.2.3)  give  the  expressions  for  Ej  {N}  and  Ex  (N).  Now 
using  the  double  exponential  functions, 


_  ....  aA'(B'-l)  +  b(B'-A') 

El{N1’  YTk 

aA,(B/  -  1)  +  b(B/  -  A') 


-i 


B'-A' 


4,4,  ,  Pji<hi 

IZPjiqiilog-1— 
Lf11=I  PjoQio 


.  PiiQn  ,  .  Puq2i 

Pu  qu  log— —  +  Pi  1  <bi  log- 


PioQio 


Pl0<l20 


.  Pnq3i  .  Pnq4i  ,  P2iQu 

+  PiiOsilog— — -  +  Pnq4ilog- - +  P2tqij1og - 

Pioqao  Pioq4o  P2oqio 


,  _  _  .  P2iq2i  .  .  P2iq3i  .  .  P2iq4i 

+  P2iq2llOg— —  +  P2iq3ilOg— —  +  P2iq4llOg— - 
P20q20  P20q30  P20q40 


.  P3iqn  .  P3iq2i  .  P3iq3i 

+  p31qulog— —  +  Psiq^ilog— — —  +  P3iq3il°gT~; — 

Psoqio  P30q20  P30q30 


.  _  _  p3iq4i  .  .  p4iqu  ,  p4 1 Q21 

+  P3iq4i log +  p4iqnlog_  - .  +  p4iq2ilog- 


p3oq4o 


p-wqio 


P40q20 


.  P4iq3i  .  P4iq4i 

+  P4iq3il°g~ — ~  +p4iq4ilog' 


P40q30 

aA'CB"  -  1)  +  b(B'  -  A') 
B'-A' 


p4oq4o 


-1 
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Ti  -Di 


Ti  -Di  .  Tj-Ih  . 

J  ±Le->-. ■  j  iic-‘-i'™,log_ri 


J  j  ~-e 


T1-D2 

l:2 


1  -X11  X-IH21  f  ^1  .-X]  I  x-jil2l 


L  2 


Ti-D,  T2-Di 

J  J  *2  c~X2!  y-tl2ll 


Ti  -Di  -  T2  . 

■  J  h***"-"*'  J 

—  1  T2-D2  Z 


Ti  -Dj  .  T2  +  D2  . 

f  Xi  I  *— H12I  (  *2  -^ly-tigl. 


L  2 


T,-D, 


^■1 

J  2 


-Xi  1  x-m2i 


It 


^2  .-ii'y- 1122 1, 


Ti-D, 


T2  +  Dl 


T2-Di, 

1  f  -± 

J  2 


T>  1  T2 

2 


f  A' 

+  J 


T, -D, 


T2-D2 


log- 


Tl  y  Ti  +  152 

I  f  i.' 


T,-Dj 


log- 


Ti  -  Di 

^  — i-i  1  x  — M-12 

t2 

"I 

[ 

1  J 

_2_  -X2I  y-Ral 

J 

2 

T2-Di 

2 

Ti-D, 

t2 

f 

A>  -X-i  1  x-|lu 

1  J 

M  -Xilx-mil 

J 

2 

T2-Di 

2 

T,-D, 

f 

^•1  -X,  1  x-(i12 

t2  +  d2 

!  f 

-X2ly-H22l 

J 

2 

t2 

2 

Ti-D! 

f 

^  -Xi  1  X-H11 

t2  +  Di 

1  r 

*■2  _-X2l  y-|i2i  1 

2  6 

J 

t2 

2  6 

Ti-D, 

f 

-Xi  1  x  — (Xl2 

1  jf 

^2  -X2ly-|122l 

J 

2 

T2+D2 

2 

Ti  -  Di 

f 

c~i|  1  *-Hu 

1  J 

^2  -X2i  y-mi 1 

J 

2 

t2+P2 

2 

Ti 

j 

e~Xl 1 

Ti-D2 

1  f 

^2  -il>  y-HH1 

T,-D, 

2 

J 

2 

Ti 

J 

^  -ii  1  x-(ln 

Ti 
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(4.4.7) 


Combining  equations  (4.4.6)  and  (4.4.7),  G  =  Ej  (N)  +  Eq{N)  can  be  calculated. 
Below  is  the  list  of  the  probability  of  each  quantized  level  under  different  Tj,  Dj,  T2, 
and  D2  where  M-n  <  Tt  <  p.12  and  |i21  <  T2  <  |i22.  It  has  been  found  by  our  evalua¬ 
tion  that  G  has  its  minima  at  the  region  jin  <  Tj  <  [ii2  and  |i2i  <  T2  <  1122. 
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J  ^Te  -  Te  -  ye 

5-D2  22  2 


Q20  =  J 

T2-D2  2 


T2_D2  2 

ql0  =  J  ^.e-^2iy-H2ii  =  !  _  l^MTj-Dz-w,) 
OQ  2  2 

If  T2  +  D2  £  P22  an(i  T2  -  D2  <  P-21 
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_  r  ^2  -X-2 1  y-^22  *  _  1  -^2(T2  +  D2-R22) 

q4i  -  J  t  -j 

T2+D2  z  z 


031 


T2  +  D2  . 

_  f  ^2  — X2  1  y-H22l  _  1 _ L«~^2(U22-T2) _ 1_  ~>-2(T2  +  D2 - R22> 

-  J  26  -1  2e  "2e 


t2  , 

f  A2  -X2 1  y-|A22 1  _  1  -^2(^22-T2)  1  • -*-2(H22-T2  +  D2) 

<i2i=  ye  - T  “  7 

T2-D2  z  z  z 


qn 


T2_D2i 

_  r  4,2  c~X2  t  y  —  M-22  I  _  j  c~*-2(^22~T2  +  D2) 

J  2  2 


j-  ^-2  -X2I  y-H21  I  _  1  -^2(T2  +  D2-ti2l) 

Q40  ~  I  “T*e  “  Te 

T2  +  D2  Z  Z 

T2  +  D2, 

_  _  f  A2  — X2 1  y-H21  1  _  1  -■ -*2(T2-H2l)  1  -2-2(T2  +  D2-^2l) 

qjo-  ^  — e  -ye  ~J' 

t2  . 

_  f  A2  -X.2ly-fl2I  I  _  1  J_  -*2(M2J -T2  +  D2)  1  -^2(T2-H22) 

q20  =  I  -r~c  -  1  -  — e  -  — e 

T2-D2  Z  Z  Z 

T2-°2, 

r  k2  -X2I  y-ji2i  I  _  1  ■ -X2(|A21  -T2  +  D2) 

qio  -  J  ~2~e  -  —t 


If  T2  +  D2  <  |l22  ^2  —  ^2  —  M-21 


_  f  ^2  -X2I  y-uzzl  _  ,  J_  ->-2(^22 -T2-D2) 

q4i  -  1  ~z~e  - 1  ve 

T2  +  L>2  Z  Z 


T2  +  D2 


q3i  - 


i 


^2  x.2  I  y — M-22  I  _  J_  -*-2<H22 -T2-D2) _ Lp-^2(H22-T2) 


2 

T2 


2  6 


_  _  f  ^2  -X2  I  y  — M-22  1  _  1  ~*-2<H22-T2)  1  ~*-2( H22 -T2  +  £>2) 

q2i-_J_— e  -ye  - -t 


t2-i>2 

T2-D2 


_  f  ^2  -X-2  1  y-M-^  I  _  J_  -^2(M2-T2  +  D2) 

qn  -  J  ~e  -  ~e 
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_  [  ^2  -X2I  y-H21  I  _  J_--*-2(T2  +  D2-«l) 

Q40  “  I  ~Z~Q  -  Te 

T2  +  D2  2  2 

T2  +  D2  , 

C  K1  -X2I  y-H21  I  _  1  -*-2(T2~H2l)  1  -^2(T2  +  D2-H2l) 

qjo-  ^  — e  -Te  -Je 

t2  , 

_  (  ^2  -X2ly-H21  1  _  _L»-X2(T2-D2-M2l)  _  l-*-2(T2-H2l) 

q20  -  J  -r-e  -  ye  -  e 

T2-D2  z  2  2 

t2  —  °2  2 

(•  k2  -X2I  y-H21 1  _  1  J_~-^2(T2-D2-^2l) 

qio  -  j  — e  - 1  ~  2 e 

T2  +  I>2  2 

r  K2  -x2l  y-M-21  I  _  1  ■ -^2(T2-H2l)  j_  -^2(T2  +  D2-Wl) 

q3o-  ^  — e  -Te  -  2* 

f  -X2!y-H21  I  _  J.  -^2(T2-D2-Wl)  j  -*-2(T2-H2l) 

q20  “  J  "T"e  “  0  e  “  0  e 

T2-D2  2  2  2 


T2  -D2  . 

_  _  f  ^2  -X2ly-tt2i  I  _  ,  1  -X2(T2-D2-H2l) 

qi°-  J  ye  ->y' 


If  T2  +  D2  <  (J-22  and  T2  —  D2  <  (i2i 


_  f  ^2_e-kly-^22l  _  j  _  J_e-X2(H22-T2-D2) 
_  J  _  9  9 


T2  +  D2 


_  _  f  pi  _  1  --MM2-T2-D2)  1  ' ~^-2(>i22-T2) 

q31  _  J  e  e 


t2  +  D2 


q2i  =  j 


^2  -X2 1  y-H22  I  _  _L_-^2(H22-T2)  _  J_  -MH22-T2  +  D2) 


T2-D2 
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t2-°2, 

f  ^2  -x2l  y— |X22  I  _  1  -*2(H22-T2+D2) 

q»  -  J  Te  -  9 e 


^2  c~X2  I  y-M-21  I  _  ^.e-^-2(T2  +  D2-H2l) 

2  2 


n  -  f  2  -X2 1  y-H21 1  -  1  r 

•  q40  -  I  -r-e  -  — e 

T2  +  D2  2  2 
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T2  +  D2}  .  . 

Kl  c~A-2  1  y-M-21  I  _  J_e-^-2(T2-Wl)  _  i_e~>-2(T2  +  D2-(i2l) 

2  2  2 

t2  1 

_  f  K1  -X2I  y-Wl  I  _  1  1  -^2(M1  -T2  +  D2)  1  -^2(T2-H22) 

420  -  I  “r~e  -  1  “  *re  -  —  e 

T2-D2  z  z  z 

T2-D2, 

_-*2ly-H21  1  _  JL(,-^2(|i21-T2  +  D2) 

_  c  —  _  c 

2  2 

4.4.  Numerical  Evaluation  of  the  G  Function 

The  numerical  evaluations  are  done  by  using  the  G  function  defined  in  the  previ¬ 
ous  section.  Since  we  have  to  deal  with  four  variables,  namely  Di ,  D2,  Ti ,  and  T2,  it 
is  not  very  convenient  to  graph  G  in  terms  of  these  variables.  The  graphs  plotted  here 
are  the  ones  with  optimal  Tt  and  T2  which  means  that  G  is  the  smallest  at  that  pairs 
of  Ti  and  T2.  In  Figure  4.2,  the  curves  represent  the  G  of  the  four-level  quantizer 
sequential  system  (FQSS)  with  respect  to  D*  and  D2.  It  is  under  the  conditions  that 
T1=T2  =  0.0,  p-n  =  -1,  ji2i  =  1.  M-i2  =  -2,  H22  =  2,  and  a  =  (3  =  0.01.  In  the  evalua¬ 
tion,  G  is  also  calculated  using  different  pairs  of  Ti  and  T2,  and  it  is  shown  in  Figure 
4.3. 

In  Figure  4.4,  the  means  are  kept  the  same  as  Figure  4.2,  and  both  a  and  P 
increased  to  0.05  with  Tj  =0.0  and  T2  =0.0.  Next,  it  shows  G  function  with  dif¬ 
ferent  error  probabilities,  a  =  0.01  and  P  =  0.05  with  Tt  =T2=0.0  in  Figure  4.5. 
With  means  shifted  to  the  left,  jin  =-2,  (i2i  =  1.  H12  =  -3,  and  ji22  =2,  the  G  is 
drawn  in  Figure  4.6  with  the  same  error  probabilities  as  the  previouse  graph  and  Tj  = 
T2  =  -0.5. 


qio  =  J 


<130  =  j 
*2 
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Figure  4.7  shows  what  happens  to  G,  when  the  means  arc  fin  =-3,  =  3, 

JJ.J2  =  -5,  and  ja.22  =  5.  The  error  probabilities  a  =  0.01  and  (3  =  0.05  with  Tj  =  -0.1 
and  T2  =  -0.2. 
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D!  =0.1 
D,  =  2.0 


Dj  =  1.5 


D,  =  0.3 


D!  =  1.0 


D,  =  0.5 
D,  =  0.7 


5 


D,  =0.8  for  T,  =  0.3  «t 


D,  =  0.7  for  T,  =-0.2  & 
D,  =0.7  fort,  =0.0* 
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Figure  4.7  G  of  FQSS  with  T1=-.lT2=-.2  and  means  at  -3,3  and  -5,5 
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CHAPTER  5 


Comparison  of  the  Performance  of  Sequential  Systems 


5.1.  Performance  of  Individual  System 

G  function  has  been  defined  as  the  sum  of  the  expected  number  of  data  to  be 
taken  by  the  system  under  hypotheses  Hi  and  Hq.  By  minimizing  G,  it  also  minim¬ 
izes  the  data  that  satisfies  both  hypotheses.  Our  optimal  sequential  system  is  obtained 
by  finding  the  minimum  of  that  function. 

In  the  previous  chapters,  the  G  function  is  evaluated  over  different  quantizers, 
means,  and  error  probabilities.  Here,  those  functions  are  compared  with  each  other 
such  that  different  characteristics  of  the  sequential  systems  are  studied. 

The  plots  that  we  discuss  here  are  from  Chapter  3  and  4. 

5.1.1.  G  of  Two-level  Quantizer  System 

In  Chapter  3,  Figure  3.2  shows  the  G  function  with  respect  to  Tj  and  T2.  As  Tj 
and  T2  come  closer  to  the  origin,  G  decreases.  It  has  its  minima  at  Tj  =0.0  and 
T2  =0.0,  and  G  =  45.26.  In  this  case,  it  seems  that  Tj  and  T2  are  independent  of 
each  other  that  no  matter  what  Tj  is,  the  minima  always  occurs  at  T2  =  0.0.  Figure 
3.3  shows  that  G  decreases  as  the  error  probabilities  increase.  This  can  be  easily 
explained  since  the  system  allows  a  larger  percentage  of  error,  it  requires  less  data  to 
be  collected  before  making  a  decision.  Again  when  Tj  and  T2  are  at  the  origin,  G 
has  its  minima  which  is  equal  to  26.63.  Figure  3.4  uses  different  error  probabilities, 
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and  a  =  0.01  and  (5  =  0.05.  The  minimum  of  G  is  equal  to  35.61  at  T;  =  T2  =  0.0. 


.  Table  5.1  lists  the  the  minimum  of  G  from  Figure  3.2,  Figure  3.3,  and  Figure 
3.4.  It  is  clear  that  G  decreases  as  the  error  probabilities  increase. 


A  little  more  interesting  observation  can  be  found  when  the  means  move  away 
from  the  origin  as  shown  in  Figure  3.6  and  3.7.  The  minimum  of  G  does  not  occur  at 
Ti  =  T2  =  0.0.  Instead,  as  in  Table  5.2,  it  has  its  minima  at  Tj  =  -0.1  and  T2  =  -0.2, 
and  a  >  [3.  The  data  of  Table  5.4  are  from  Figure  3.6.  From  the  table,  one  can  notice 
that  T2  at  the  minimum  G  increases  as  Ti  decreases.  Figure  3.7  shows  the  case  in 
which  the  error  probability  a  is  less  than  {3,  both  Ti  and  T2  are  shift  to  the  left,  and 
Tj  =0.3  and  T2  =  0.2.  The  minimum  value  of  G  is  6.509  which  is  smaller  than  the 
one  in  Figure  3.4  under  the  same  error  probabilities.  This  shows  that  it  is  easier  to 
detect  a  signal  when  the  two  density  functions  under  hypothesis  Hq  and  Hj  arc 
further  apart. 


Thresholds 

a=  (3  =  0.01 

a  =  0.01  (3  =  0.05 

a  =  (3  =  0.05 

T, 

t2 

G 

G 

G 

-0.1 

0.0 

49.68 

39.18 

29.23 

-0.8 

0.0 

48.02 

37.86 

28.26 

-0.6 

0.0 

46.31 

36.87 

27.53 

-0.4 

0.0 

45.93 

36.18 

27.03 

-0.2 

0.0 

45.42 

35.76 

26.73 

-0.1 

0.0 

45.30 

35.66 

26.66 

* 

O 

O 

*0.0 

45.26 

35.61 

26.63 

0.1 

0.0 

45.26 

35.64 

26.66 

0.2 

0.0 

45.42 

35.72 

26.73 

0.4 

0.0 

45.93 

36.10 

27.03 

0.6 

0.0 

46.78 

36.76 

27.53 

0.8 

0.0 

48.02 

37.72 

28.26 

1.0 

0.0 

49.68 

39.00 

29.23 

Table  5.1  Tabulated  Data  of  Minimum  G  with  Two-level  Quantizer 
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a  = 

0.05  3  = 

0.01 

T, 

T2 

G 

-0.3 

-0.1 

8.08 

-2.5 

-0.1 

7.51 

-2.0 

-0.1 

7.16 

-1.5 

-0.1 

6.82 

-1.0 

-0.1 

6.64 

-0.8 

-0.1 

6.59 

-0.6 

-0.1 

6.55 

-0.4 

-0.2 

6.53 

-0.2 

-0.2 

6.517 

*-0.1 

*-0.2 

6.516 

0.0 

-0.2 

6.519 

0.1 

-0.2 

6.52 

0.2 

-0.2 

6.53 

0.4 

-0.2 

6.56 

0.8 

-0.2 

6.64 

1.0 

-0.2 

6.71 

1.5 

-0.3 

6.93 

2.0 

-0.3 

7.24 

2.5 

-0.3 

7.68 

2.9 

-0.3 

8.14 

Table  5.2  Minimum  G  with  Two-level  Quantizer  a  =  0.05  (3  =  0.01 
5.1.2.  G  of  the  Four-level  Quantizer  System 

In  general,  as  it  shows  in  the  graphs  of  Chapter  4,  both  D1  and  D2  do  not  change 
for  a  given  means.  However,  like  previous  section,  the  G  function  decreases  as  the 
error  probabilities  increase.  In  Table  5.3,  the  information  in  Figure  4.2,  Figure  4.4, 
and  Figure  4.5  is  tabulated. 


As  in  Figure  4.7,  it  has  the  smallest  minima  of  G  compared  to  other  figures 
because  the  means  are  further  apart  from  each  other  which  makes  the  target  more  dis¬ 
tinguishable.  This  allows  the  system  to  require  less  data  to  make  the  decision 
whether  the  target  is  present.  The  data  from  Figure  4.7  is  tabulated  in  Table  5.4. 


6  A 


H 

11 

H 

Ki 

II 

O 

O 

a  =  (3  =  0.01  |  a  =  0.01  0  =  0.05 

D1 

D2  | 

a  =  0  =  0.05 

G 

G 

0.1 

1.3 

43.52 

34.24 

25.61 

0.3 

1.3 

43.00 

33.82 

25.29 

0.5 

1.3 

42.69 

33.59 

25.12 

0.6 

1.3 

42.64 

33.55 

25.09 

*0.7 

*1.3 

42.63 

33.548 

25.088 

0.8 

1.3 

42.67 

33.58 

25.11 

1.0 

1.3 

42.86 

33.73 

25.22 

1.5 

1.3 

43.27 

34.05 

25.46 

2.0 

1.3 

43.49 

34.20 

25.58 

3.0 

1.3 

43.65 

34.45 

25.69 

4.0 

1.3 

43.74 

34.42 

25.74 

Table  5.3  Tabulated  Data  of  Minimum  G  with  Four- level  Quantizer 


T,  = -0.1  T2  = -0.2 

a  =  0.05  0  =  0.01 

D1 

D2 

G 

0.1 

3.3 

6.14 

0.5 

3.3 

6.78 

1.0 

3.3 

5.84 

1.5 

3.3 

5.76 

1.7 

3.3 

5.75 

1.9 

3.3 

5.743 

*2.0 

*3.3 

5.742 

2.1 

3.3 

5.744 

2.3 

3.3 

5.75 

2.5 

3.3 

5.77 

3.0 

3.3 

5.84 

4.0 

3.3 

5.97 

5.0 

3.3 

6.04 

Table  5.4  Minimum  G  with  Four-level  Quantizer,  a  =  0.05  p  =  0.01 
Table  5.5  tabulates  the  data  of  various  Tj  and  T2  with  a  =  0.01,  p  =  0.01, 
un  =  —1,  U12  =  1,  U21  =  -2,  and  U22  =  2.  This  shows  how  G  varies  with  T!  and  T2. 


5.2.  Comparison  of  System 

From  the  numerical  analysis,  one  can  conclude  that  the  system  with  a  four-level 
quantizer,  in  average,  requires  less  data  to  make  a  decision  than  the  one  with  a  two- 

level  quantizer.  By  compared  by  picking  a  column  from  both  Table  5.1  and  Table 
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T,  =  -0.2  T2  =  -0.3 

o 

o 

ii 

It 

H 

T,  =  0.3  T2  =  0.5 

D2 

D1 

G 

D1 

G 

D1 

G 

0.1 

0.7 

43.84 

0.7 

43.70 

0.8 

44.06 

0.3 

0.7 

43.49 

0.7 

43.36 

0.8 

43.69 

0.5 

0.7 

43.21 

0.7 

43.09 

0.8 

43.39 

0.6 

0.7 

43.10 

0.7 

42.98 

0.8 

43.27 

0.8 

0.7 

42.92 

0.7 

42.81 

0.8 

43.08 

1.0 

0.7 

42.80 

0.7 

42.70 

0.8 

42.94 

1.5 

0.7 

42.73 

0.7 

42.65 

0.8 

42.84 

2.0 

0.7 

42.90 

0.7 

42.87 

0.8 

43.94 

3.0 

0.7 

43.33 

0.7 

43.27 

0.8 

43.39 

4.0 

0.7 

43.55 

0.7 

43.48 

0.8 

43.67 

Table  5.5  G  with  Various  Tj  and  T2 

5.3.  G  of  the  two-level  system  is  larger  than  the  one  of  the  four-level  system.  For 
a  =  0.01,  p  =  0.01,  uii  =  -1,  u12  =  1,  u2j  =  -2,  and  u22  =  2,  the  improvement  on  the 
G  function  is  5.8  %.  However,  for  a  =  0.05,  P  =  0.01,  un  =-3,  ui2  =  3,  u2J  =  -5, 
and  u22  =  5,  the  improvement  on  the  G  function  is  11.88  %.  Table  5.6  shows  the 
improvement  on  the  G  function  for  various  error  probabilities  and  means. 

Notice  from  Table  5.6,  the  improvement  of  G  for  a  given  set  of  means  are  the 
same  no  matter  what  the  error  probabilities. 

The  other  observation  is  that  Ti  and  T2  of  the  four-level  quantizer  system  is  the 
same  of  the  two-level  quantizer  system.  Ti  andT2  from  a  two-level  quantizer  system 
can  be  used  for  designing  a  four-level  quantizer  if  it  is  under  same  means  and  error 
probabilities.  This  approach  will  reduce  the  total  number  of  calculation  to  find  the 
optimal  G. 


Condition 

G  with 

Two-level  Quantizer 

G  with 

Four-level  Quantizer 

%  Improvement 

a  =  0.01  3  =  0.01 
tin  =  — 1  u,2  =  1 

u2i  =  -2  U22  =  2 

45.26 

42.63 

5.8% 

a  =  0.01  (3  =  0.05 
Uu  =-l  u12=  1 

U21  =  —2  U22  =  2 

26.63 

25.09 

5.8% 

a  =  0.05  P  =  0.05 
u  1 1  =-l  Ui2=  1 

U21  =  —2  U22  =  2 

35.61 

33.55 

5.8% 

a  =  0.05  P  =  0.01 
uu  =  -3  u12  =  3 
u2i  =  -5  u22  =  5 

6.516 

5.742 

n.88% 

Table  5.6  Improvement  on  G 
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CHAPTER  6 


System  Simulations 


6.1.  Simulation  Method 

From  Chapter  3  and  4,  the  optimal  sequential  detection  systems  are  found.  It 
has  been  shown  that  the  quantized  levels  of  a  quantizer  has  a  major  effect  on  the 
sequential  system.  By  selecting  the  quantizers  that  associate  with  the  optimal  system, 
the  system  can  be  evaluated  under  the  real  environment. 

6.1.1.  Derivation  of  The  Double  Exponential  Environment 

Since  the  computer  that  is  used  to  simulate  the  system  has  a  random  value  gen¬ 
erator  with  uniform  distribution,  transformation  method  must  be  used  to  generate  the 
double  exponential  density  environment. 

Let  x  be  the  uniform  distributed  random  variable,  and  x  e  (0,  1),  then  the  double 
exponential  density  function  Fy(y)  is  equal  to  x.  Thus  y  =  Fy-1  (x). 

The  distribution  function  can  be  evaluated  by  integration  of  the  density  function 
fy(x)  =  ye"Xlx_^l,and 

Fy(x)=  Jf  ye’xiz~^dz 

If  x  >  (i,  then 

Fy(x)=  1  - 


68 


If  x  <  (j.,  then 


the  inverse  fuction  Fy(y)  is 


Fy(y)  =  ye 


y  =  |i  +  ~log(  2x)  for  x  >  0.5 

A, 

y  =  M-  —  ^-l°g(  2(1  -x))  for  x  <  0.5 


6.1.2.  Simulation  and  Discussion 

The  first  systems  to  be  simulated  are  under  the  conditions  that  p.n  =  -1,  p12  =  1, 
(j-2i  =  -2,  and  |i.22  =  2.  With  the  error  probabilities  varing,  the  amount  of  data  N 
required  to  make  a  decision  are  collected.  Total  of  one  thousand  trials  are  made,  and 
Figure  6.1  shows  N  of  two-level  quantizer  system  with  error  probabilities  a  =  0.01 
and  (3  =  0.01.  Figure  6.2  is  N  of  four- level  quantizer  system  with  same  error  probabil¬ 
ities  as  in  Figure  62.  In  Figure  6.3,  Figure  6 .4,  and  Figure  6  .5,  it  gives  an  enlarged 
view  on  how  N  varies  in  different  trial  and  error  probabilities. 


Figure  6.6  assumes  that  the  means  |Xn  =  ~3,  ji12  =  3,  jx21  =  -5,  and  |i22  =  5,  and 
N  under  both  quantized  systems  are  plotted.  Notice  that  the  dashed  line  of  all  plots 
represents  N  of  the  four-level  quantizer  system. 

From  the  graphs,  one  may  notice  that  the  system  with  four-level  quantizer  has  a 
higher  G  values  than  the  system  with  two-level  one.  The  results  of  the  system  simu¬ 
lations  confirm  that  in  average  the  four-level  quantizer  system  requirs  less  amount  of 
data  than  the  two-level  one. 
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Trial  Number 

Figure  6.4  Enlarged  N  with  a  =  0.01  {3  =  0.05 
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CHAPTER  7 


Conclusion 

The  objective  of  this  study  is  to  analysis  and  design  a  two-sensor  sequential 
detection  system.  The  optimal  system  is  found  by  numerical  analysis,  and  that  sys¬ 
tem  is  simulated. 

In  this  study,  the  expected  number  of  observation  is  first  derived,  and  then  func¬ 
tion  G  is  defined  as  the  sum  of  that  expected  number  under  both  hypotheses.  By 
minimizing  G,  it  minimizes  the  expected  value  of  number  of  observations  required. 

G  is  directly  related  to  the  quantizer  of  each  sensor.  As  shown  in  chapter  2, 
function  G  is  an  equation  in  terms  of  the  probability  of  the  quantized  levels. 
Minimum  G  can  be  found  by  selecting  appropriate  quantizer. 

Two-level  and  four-level  quantized  system  are  used  in  the  analysis.  Assuming 
the  environment  is  double  exponentially  distributed,  G  is  evaluated  under  various 
conditions.  The  numerical  evaluations  show  that  a  four-level  quantizer  system  has 
smaller  value  of  G  than  a  two-level  one.  This  shows  that  a  four-level  quantizer  sys¬ 
tem  requires  less  data  than  a  two-level  one  to  make  a  decision. 

It  is  also  found  that  for  a  given  density  function,  the  quantizers  are  always  the 
same  for  both  sensors  to  have  a  minimum  value  of  G  no  matter  what  are  the  error 
probabilities.  As  the  means  of  the  density  function  moving  further  apart,  G 
decreases,  and  the  quantized  levels  are  not  necessarily  divide  the  probability  space 
symmetrically.  The  system  simulation  confirmed  the  results  of  the  evaluation. 


In  the  future,  it  might  be  interesting  to  look  into  the  case  of  dependent  observa¬ 


tion. 
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APPENDIX  A 


Program  for  Evaluation  of  G  Function  of  Two-level  Quantizer  System 


This  program  calculates  the  G  function  of  two-level  quantizer  system  for  vari¬ 
ous  parameters.  The  means  of  double  exponential  and  error  probabilities  are  pro¬ 
vided,  and  then  this  program  calculate  the  probabilities  associated  with  each  quan¬ 
tized  level  for  different  two-level  quantizers.  The  main  routine  find  the  value  of  G. 

The  plots  in  chapter  3  are  provided  by  this  program. 
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This  program  is  used  to  find  the  minimun  number  of  data  that 
required  to  determine  wheter  a  target  is  present  by 
appling  sequencial  detection  theory.  Two  independent  sensers 
are  used  to  detect  the  signals  from  the  source.  Two-level 

quantizers  are  used  here,  and  optimal  quantizer  is  found. 
********•*****«**•****«**»•*•*•**•••*••*•********••*•••»*******••••*/  #include 

<stdio.h>  #include  <math.h>  /"define  the  value  of  error  probabilities*/  #define  alpha  0.05  #define  beta 
0.01  /"define  the  constants  of  double  exponential  density  functions*/  #define  U10  -3  #define  Ull  3 
#dcfine  U20  -5  #define  U21  5  /define  11  0.25  #define  12  0.125 


mainO  /*  find  min  EN  and  store  the  data  for  further  uses*/  { 
double  T1  ,T2;  /"quantized  levels  of  both  sensors*/ 
double  ConslO.ConsOO; 
double  U10.U20; 
double  P10.P20: 
double  f; 
double  m,ta,tb; 
double  U0; 

FILE  *fpout; 

fpout  * 


* 


fopen("chp3fg6"  t"w"); 
calculate  G  func¬ 


tion 


*  ************»************************************»***********^ 


for  ( T1  =  -3.0;  T1  <=  3.0;  T1  =  T1  +  0.1)  { 
for  (  T2  =  -5.0;  T2  <=  5.0;  T2  =  T2  +  0.1)  { 
f  =  (ConslO  /U1(T1,T2))  +  (ConsOO  /  U2(T1,T2)); 
fprintf(fpout,"%5.2f  %7.2f  %10.4fD,Tl,T2L); 

} 


) 

fclose(fpout); } 

double  U(pl,pO,ql,qO,cl,c2)  /*  calculate  the  U  expression  */ double  pl,p0,ql,q0;  double  cl,c2;  { 
double  i; 

i=  cl  *log((pl*(l-p0))/(p0*(l  -pi))) 

+  c2  *  log((ql  *  (1  -q0))  /  (qO  *  (1  -  ql))) 

+  log(((l  -  pi)  *  (1  -  ql))  /  ((1  -  p0)  *  (1  -  q0))); 
retum(i);  } 

/*  find  the  probabilities  of  quantized  levels*/  double  Pl(l,t,u)  double  l,t.u;  { 
double  i; 

i  =  0.5  *  exp((-l)  *  (t  -  u)); 
retum(i);  ) 


double  P2(l,t,u)  double  l,t,u;  { 


double  i; 

i  =  1  -  0.5  *  exp(  (-1)  *  (u  - 1)); 
retum(i); ) 

r  find  the  value  of  U1  which  involves  prob  p  &  q  */  double  Ul(tl,t2)  double  tl,t2;  { 
double  pl,p0,ql,q0; 
double  cl,c2; 
double  P10J*20: 
double  U0; 

double  ull,ul0,u21,u20; 

double  LI  JL2; 

ull  =U11; 

ulO  =  U10; 

u2l  =  U21; 

u20  =  U20; 

Li  -11; 

L2  =12; 

if  (tl  >=  ul  1  &&  t2  >=  u21)  { 
pi  =Pl(Ll,tl,ull); 
pO  =  Pl(Ll,tl,ulO); 
ql  =  Pl(L2,t2,u21); 
qO  =  Pl(L2,t2,u20); 
cl  =  pi; 
c2  =  ql; 

return(  U(pl  ,p0,ql  ,q0,c  1  ,c2)); 

} 

if  (tl  >=  ul  1  &&  t2  >=  u20  &&  t2  <  u21)  { 
pl=Pl(Ll,tl.ull); 
p0  =  Pl(Ll,tl,ulO); 
ql  =P2(L2,t2,u21); 
qO  =  Pl(L2,t2,u20); 
cl  =  pi; 
c2  =  ql; 

rctum(  U(p  1  ,p0,ql  ,q0,c  1  ,c2)); 

} 

if  (tl  >=  ull  &&  t2  <  u20)  { 
pi  =Pl(LUl.ull); 
pO  =  PI(LI,tI,ulO); 
ql  =  P2(L2,t2,u21); 
qO  =  P2(L2,t2,u20); 
cl  =  pi; 
c2  =  ql; 

retum(  U(p  1  ,pO,q  l  .qO,c  1  ,c2)); 

} 
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if  (tl  <  ull  &&  tl  >=  ulO  &&  t2  >=  u21)  { 
pi  =P2(LUl.ull); 
pO  =  Pl(LUl,ulO); 
ql  =  P1(L242,u21); 
qO  -  Pl(L2j2,u20); 
cl  =pl; 
c2  =  ql; 

ietum(  U(pl4i04ll^0^1^2)); 

} 

if  (tl  <  ull  &&  tl  >=  ulO  &&  t2  <  u21  &&  t2  >=  u20)  { 
pi  =P2(Ll4l.ull); 
pO  =  Pl(LUJ.ulO); 
ql  =P2(L242,u21); 
qO  =  P1(L2o2,u20); 
cl  =  pi; 
c2  =  ql; 

ieturn(  U(p  1 4)0^1  14)04:  14:2)); 

} 

if  (tl  <  ull  &&ll  >=  ulO&A  t2  <  u20)  { 
pi  =P2(LUl,ull); 
pO  =  Pl(LUl,ulO); 
ql  -  P2(L2j2,u21); 
qO  =  P2(L242,u20); 
cl -pi; 
c2  =  ql; 

retum(  U(pl  4)0 4)1 4)04:  1 4:2)); 

} 

if  (tl  <  ulO  &&  t2  >=  u21)  { 
pi  =P2(LUl,ull); 
pO  =  P2(Ll4l,ulO); 
ql  =Pl(L242.u21); 
q0  =  Pl(L242,u20); 
cl  -pi; 
c2  =  ql; 

reuim(  U(p  1 4)041  l4}0x  1^2)); 

J 

if  (tl  <  ulO  <fe&  t2  <  u21  &&.  t2  >=  u20)  { 
pi  =P2(Ll4l,ull); 
pO*P2(Ll4l.ulO); 
ql  *  P2(L242,u21); 

'  q0  =  Pl(L242,u20); 
cl  «pl; 
c2-ql; 

ietum(U(pl4)04ll  4)04:1 4:2)); 

} 
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if  (tl  <  ulO  &&  t2  <  u20)  { 
pi  =P2(LUl,ull); 
pO  =  P2(LUl,ulO); 
ql  =  P2(L2,t2,u21); 
qO  =  P2(L2,t2,u20); 
cl  =pl; 
c2  =  ql; 

retum(  U(pl,pO,ql  ,qO,c  1  ,c2)); 

}} 

I*  find  the  value  of  U2  in  function  G  which  involves  prob  p , 
double  pl,pO,ql,qO; 
double  cl, c2; 
double  P10.P20; 
double  UO; 

double  u  1 1  ,u  10,u2 1  ,u20; 

double  LI  JL2; 

ull  =  U11; 

ul0  =  U10; 

u21  =U21; 

u20  =  U20, 

LI  =  11; 

U  =12; 


if  (tl  >=  ul  1  &&  t2  >=  u2l)  { 
pi  =Pl(Ll,tl,ull); 
pO  =  Pl(Ll,tl,ulO); 
ql  =Pl(L2,t2,u21); 
qO  =  Pl(L2,t2,u20); 
cl  =p0; 
c2  =  qO; 

retum(  U(p  1  ,pO,q  1  ,qO,c  1  ,c2)); 

) 

if  (tl  >=  ul  1  &&  t2  >=  u20  &&  t2  <  u21)  { 
pi  =Pl(Ll,tl,ull); 
pO  =  Pl(Ll,tl,ulO); 
ql  =P2(L2,l2,u21); 
qO  =  Pl(L2,t2,u20); 
cl  =p0; 
c2  =  qO; 

retum(  U(p  1  ,pO,q  1  ,qO,c  1  ,c2)); 

} 

if  (tl  >=ull  &.&.  t2  <  u20)  { 
pi  =  Pl(Ll,tl,ull); 
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q  */  double  U2(tl,t2)  double  tl,t2;  { 


pO  =  Pl(Ll,tl,ulO); 
ql  =P2(L2,t2,u21); 
qO  =  P2(L2,t2,u20); 
cl  =pO; 
c2  =  qO; 

retum(  U(pl,p0,ql,q0,cl,c2)); 

} 

if  (tl  <  ull  &&  tl  >=  ulO  &&  a  >=  u21)  ( 
pi  =P2(Ll,tl,ull); 
pO  =  Pl(Ll,tl,ulO); 
ql  =Pl(L2,l2,u21); 
qO  =  PI(L2,t2,u20); 
cl  =pO; 
c2  =  qO; 

retum(  U(pl,p0,ql,q0,cl,c2)); 

} 

if  (tl  <  ull  &&tl>=ul0&&t2<u21  &&  t2  >=  u20)  { 
pi  =P2(Ll,tl,ull); 
pO  =  Pl(Ll,tl,ulO); 
ql  =  P2(L2,t2,u21); 
qO  =  Pl(L2,t2,u20); 
cl  =pO; 
c2  =  qO; 

retum(  U(pl,p0,ql,q0,cl,c2)); 

} 

if  (tl  <  ull  &&  tl  >=  ulO  &&  t2  <  u20)  { 
pi  =P2(Ll,tl,ull); 
pO  =  Pl(Ll,tl,ulO); 
ql  =  P2(L2,t2,u21); 
qO  =  P2(L2,t2,u20); 
cl  =pO; 
c2  =  qO; 

rctum(  U(pl,pO,ql,qO,cl,c2)); 

} 

if  (tl  <  ulO  &&  12  >=  u21)  ( 
pi  =P2(Ll,tl,ull); 
pO  =  P2(Ll,tl,ulO); 
ql  =  Pl(L2,t2,u21); 
qO  =  Pl(L2,t2,u20); 
cl  =pO; 
c2  =  qO; 

retum(  U(pl,p0,ql,q0,cl»c2)); 

} 

if  (tl  <ul0&&t2<u21  &&  t2  >=  u20)  { 
pi  =P2(Ll,tl,ull); 
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pO  =  P2(Ll,tl,ulO); 
ql  =  P2(L2,t2,u21); 
qO  =  Pl(L2,t2,u20); 
cl  =pO; 
c2  =  qO; 

retum(  U(pl,p0,ql,q0,cl,c2)); 

} 

if  (tl  <  ulO  &&  l2  <  u20)  { 
pi  =  P2(Ll,tl,ull); 
pO  =  P2(Ll,tl,ulO); 
ql  =  P2(L2,t2,u21); 
qO  =  P2(L2,t2,u20); 
cl  =pO; 
c2  =  qO; 

rctum(  U(pl,p0,ql,qO.cl,c2)); 

}  ) 

/* calculate  constants  in  function  G*/  double  Cons  10  { 
double  a,b; 
double  A3; 
double  i; 

A  =  (beta)  /  (1  -  alpha); 

B  *  (1  -  beta)/  alpha; 
a  =  log(A); 
b  =  log(B); 

i-(a*A*(B-l)  +  b*B*(l-  A))  /  (B  -  A); 
retum(i); ) 

double  ConsOO  { 
double  a,b; 
double  A3; 
double  i; 

A  =  (beta)  /  (1  -  alpha); 

B  =  (1  -  beta)  /  alpha; 
a  =  log(A); 
b  =  log(B); 

i  =  (a  *  (B  -  1)  +  b  *  (1  -  A))  /  (B  -  A); 
retum(i);  ) 
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APPENDIX  B 


Program  for  Evaluation  of  G  Function  of  Four-level  Quantizer  System 


This  program  calculates  the  G  function  of  four-level  quantizer  system  for  vari¬ 
ous  parameters.  The  means  of  double  exponential  and  error  probabilities  are  pro¬ 
vided,  and  then  this  program  calculate  the  probabilities  associated  with  each  quan¬ 
tized  level  for  different  four-level  quantizers.  The  main  routine  find  the  value  of  G. 

The  plots  in  chapter  4  are  provided  by  this  program. 
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This  program  is  used  to  find  the  minimum  number  of  data  that 
required  to  determine  whether  a  target  is  present  by 
appling  sequencial  detection  theory.  Two  independent  sensers 
are  used  to  detect  the  signals  from  the  source.  G  function  is 

calculated  using  various  four-level  quantizers. 

** ** **************************************************************** j  #include 

<stdio.h>  ^include  <math.h>  #define  alpha  0.05  #define  beta  0.01  #define  U10  -3  #define  Ull  3 
tfdefine  U20  -5  #define  U21  5  #define  11  0.25  #define  12  0.125  #define  T1  -0.1  #define  T2  -0.0 

mainO  /*  find  min  G  and  store  the  data  for  further  uses*/  { 
double  D1,D2; 
double  ConslO.ConsOO; 
double  U10; 

double  P10,P20,P30,P40,P50  WO; 

double  f; 
double  m,d3,d4; 
double  UO; 
int  k; 

FILE  *fpout,*fp; 


m  =  1000.0; 
d3  =  0.0; 

d4  =  0.0;  ft********************************************************  *  calculate  G 
function  * •****«********************************•******************/ 

fpout  =  fopen("chap4fig6.t","w"); 
fp  =  fopen("supfig6.t”,"w"); 
for  ( D1  =  0.1;  D1  <=  3.0;  D1  =  D1  +0.1)  { 
for  ( D2  =  0.1;  D2  <=  5.0;  D2  =  D2  +  0.1)  { 
f  =  (Cons  10  /  U1(D1,D2.1»  +  (ConsOO  /  U1(D1,D2,0)); 
if  (  f  <=  m)  { 
m  =  f; 
d3  =  Dl; 
d4  =  D2;  . 

) 

fjprintf(fpout,"%5.2f  %7.2f  %10.4fO.Dl.D2.0; 

) 

) 

fprintf(fp,”%5.2f  %5.2f  %10.4f0,d3,d4.m); 

fclose(fpout); 

fclose(fp);  ) 


j  *■*■************+***»  ************  ********  *******  ******  *•*****»**•« 


double  U(p41,p31,p21,pll,p40,p30,p20,pl0,q41,q31,q21,qll,q40,q30,q20,ql0jc)  /*  calculate  the  U 


expression  */ double  p41,p31,p21,pll,p40,p30,p20,pl0;  double  q41,q31,q21,qll,q40,q30,q20,ql0;  int 
k;  { 

double  i,a,b; 

‘  if(k  —  0){ 

a=  plO  *  qlO  *  log((pll  *  qll)/(plO  *  qlO)) 

+  p20  *  qlO  *  log((p21  *  ql  1)  /  (p20  *  qlO)) 

+  p30  *  qlO  *  log((p31  *  ql  1)  /  (p30  *  qlO)) 

+  p40  *  qlO  *  log((p4 1  *  qll)  /  (p40  *  qlO)) 

+  plO  *  q20  *  log((pl  1  *  q21)  /  (plO  *  q20)) 

+  p20  *  q20  *  log((p21  *  q21)  /  (p20  *  q20)) 

+  p30  *  q20  *  log((p3 1  *  q21)  /  (p30  *  q20)) 

+  p40  *  q20  *  log((p41  *  q21)  /  (p40  *  q20)); 
b  =  plO  *  q30  *  log((pl  1  *  q31)  /  (plO  *  q30)) 

+  p20  *  q30  *  log((p21  *  q3 1)  /  (p20  *  q30)) 

+  p30  *  q30  *  log((p31  *  q31)  /  (p30  *  q30)) 

+  p40  *  q30  *  log((p41  *q31)/(p40*  q30)) 

+  plO  *  q40  *  log((pl  1  *  q41)  /  (plO  *  q40)) 

+  p20  *  q40  *  log((p21  *  q41)  /  (p20  *  q40)) 

+  p30  *  q40  *  log((p31  *  q41)  /  (p30  *  q40)) 

+  p40  *  q40  *  log((p41  *  q41)  /  (p40  *  q40)); 
i  =  a  +  b; 
retum(i); 

} 

if(k=l){ 

a  =  pi  1  *  ql  1  •  log((p  1 1  *  ql  1)  /  (plO  *  qlO)) 

+  p21  *  qll  *  log((p21  *  ql  1)  / (p20  *  qlO)) 

+  p31  *  ql  1  *  log((p31  *  ql  1)  /  (p30  *  qlO)) 

+  p41  *  qll  *  log((p41  *  qll)/(p40  *  qlO)) 

+  pi  1  *  q21  *  log((pl  1  *  q21)  /  (plO  *  q20)) 

+  p21  *  q21  *  log((p21  *  q21)  /  (p20  *  q20)) 

+  p31  *  q21  *  log((p31  *  q21)  /  (p30  *  q20)) 

+  p41  *  q21  *  log((p41  *  q21)  /  (p40  *  q20)); 
b  =  pi  1  *  q31  *  log((pl  1  *  q31)  /  (plO  *  q30)) 

+  p21  *  q31  *  log((p21  *  q31)  /  (p20  *  q30)) 

+  p31  *  q31  *  log((p3 1  *  q31)  /  (p30  *  q30)) 

+  p41  *  q31  *  log((p41  *  q31)  /  (p40  *  q30)) 

+  pi  1  *  q41  *  log((pl  1  *  q41)  /  (plO  *  q40)) 

+  p21  *  q41  *  log((p21  *  q41)  /  (p20  *  q40)) 

+  p31  *  q41  *  log((p3 1  *  q41)  /  (p30  *  q40)) 

+  p41  *  q41  *  log((p41  *  q41)  /  (p40  *  q40)); 
i  =  a  +  b; 
retum(i); 

}) 


/ft*****************************************************/ 
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/*  find  die  probabilities  of  quantized  levels  */  double  Pl(u,l,t)  double  l,t,u;  [ 

double  i; 

i  =  0.5  *  exp((-l)  *  (t  -  u)); 
retum(i); } 

double  P2(u,l,t,d)  double  uj,t,d;  { 
double  i; 

i  =  1  -  0.5  *  exp((-l)  *  (u  - 1))  -  0.5  *  exp((-l)  *  (t  +  d  -  u)); 
return(i);  } 

double  P3(u4,t)  double  l,t,u;  { 
double  i; 

i  *  0.5  *  exp((-l)  *  (u  - 1)); 
rctum(i);  } 

double  P4(uJ,t)  double  l,t,u;  ( 
double  i; 

i  =  1  -  0.5  *  exp(  (-1)  *  (u  - 1)); 
retum(i); ) 

double  P5(uJ,t,d)  double  uj,t,d;  { 
double  i; 

i  =  -  0.5  *  exp((-l)  *  (u  - 1))  +  0.5  *  exp((-l)  *  (u  - 1  -  d»; 
retum(i); } 

double  P6(u4,t,d)  double  u  J,t,d;  ( 
double  i; 

i  =  0.5  *  exp((-l)  *  (t  -  u))  -  0.5  *  exp((-l)  *  (t  +  d  -  u)); 
retum(i);  } 

double  P7(uJ,t)  double  l,t,u;  { 
double  i; 

i  =  1  -  0.5  *  exp(  (-1)  *  (t  -  u)); 
retum(i);  ] 


/* ♦ **** ♦ ** * ** *♦* ****** * * 


it******************************  I 


I*  find  the  value  of  U1  which  involves  prob  p  &  q  */  double  Ul(dl,d24c)  double  dl,d2;  int  k;  { 
double  pll,p21,p31,p41; 
double  pl0,p20,p30,p40; 
double  q  1 1  ,q2 1  ,q3 1  ,q4 1 ; 
double  ql0,q20,q30,q40; 
double  P104>204>304,404>50,P60.P70; 
double  U0; 

double  ull,ul0,u21,u20; 
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double  LI  X2; 
doublet  1,(2; 
tl  =  Tl; 
t2  =  T2; 
ull  =  U11; 

UlO  as  UlO, 
u21  =  U21; 
u20  =  U20; 

LI  =11; 

L2  =12; 

if  ((tl+dl)  <  ul  1  &&  (tl-dl)  >=  ulO  &&  (t2+d2)  <  u21  &&  (t2-d2)  >=  u20) 

[ 

p41  =P4(ullJ-l, tl+dl); 
p31  =  P5(ullJLl,tl,dl); 
p21  =  P5(ullJLl,tl-dI,dl); 
pll  =P3(ullJLl,tl-dl); 
p40  =  Pl(ulOJLl,tl+dl); 
p30  =  P6(ulOJLl,tltdl); 
p20  =  P6(ulOJLl,tl-dl,dl); 
plO  =  P7(ulO,Ll,tl-dl); 

q41  =  P4(u21JL2,t2+d2); 
q31=P5(u21L2,t2,d2); 
q21  =  P5(u21L2,t2-d2,d2); 
ql  1  =  P3(u21,L2,t2-d2); 
q40  =  Pl(u20JJ2,t2+d2); 
q30  =  P6(u20,L2,t2,d2); 
q20  =  P6(u20L2,t2-d2,d2); 
qlO  =  P7(u20,L2,t2-d2); 

reuim(  U(p4 1  ,p3 1  ,p21  ,pl  1  ,p40,p30,p20,p  10, q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,q  10, 

k)); 

) 

if((tl+dl)<ull  &&  (tl-dl)>=  ulO  &&  (t2+d2)<u21  &&  (t2-d2)<u20  ) 

{ 

p41  =  P4(ullJLl, tl+dl); 
p31  =  P5(ul  lX-l.tl.dl); 
p21  =  P5(ulU-l,tl-dl,dl); 
pll  =  P3(ullJLl, tl-dl); 
p40  =  Pl(ul0i-l, tl+dl); 
p30  =  P6(ul0jLI,tl,dl); 
p20  =  P6(ul0Ll,tl-dl,dl); 
plO  =  P7(ulOLl, tl-dl); 
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q41  =  P4(u2U-2,t2+d2); 
q31  =  P5(u21X2,t2,d2); 
q21  -  P5(u21J-2,t2-d2,d2); 
ql  1  =  P3(u21,L2,t2-d2); 
q40  =  Pl(u20J-2,t2+d2); 
q30  =  P6(u20JL2,t2,d2); 
q20  =  P2(u20J-2,t2-d2,d2); 
ql0  =  P3(u20JL2,t2-d2); 

rctum(  U(p4 1  ,p3 1  ,p2 1  ,p  1 1  ,p40,p30,p20,p  10, q4 1  ,q3 1  ,q2 1  ,q  1 1  tq40,q30,q20,q  1 0, 

k)); 

} 

if((tl+dl)<ull  &&  (tl-dl)>=  ulO  &&  (t2+d2)>=u21  &&  (t2-d2)>=u20) 

{ 

p41  =  P4(ullJLl,tl+dl); 
p31  =  P5(ullJ-l,tl,dl); 
p21=P5(ullXl,tl-dl.dl); 
pi  1  =  P3(ul  1  JLl,tl-dl); 
p40  =  Pl(ulOJ-l,tl+dl); 
p30  =  P6(ulO,Ll,tl,dl); 
p20  =  P6(ul0^1,tl-dl,dl); 
plO  =  P7(ulO,Ll,tl-dl); 

q41  =:Pl(u21JL2,t2+d2); 
q31=P2(u21,L2,t2,d2); 
q21  =  P5(u21JL2,t2-d2,d2); 
ql  1  =  P3(u21  JL2,t2-d2); 
q40  =  Pl(u20X2,t2+d2); 
q30  =  P6(u20X2,t2,d2); 
q20  =  P6(u20X2,t2-d2,d2); 
qlO  =  P7(u20X2,t2-d2); 

retum(  U(p4 1  ,p3 1  ,p2 1  ,p  1 1  ,p40,p30,p20,p  10,q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,q  10, 

k)); 

) 

if((tl+dl)<ull  &&  (tl-dl)>=  ulO  &&  (t2+d2)>=u21  &&  (t2-d2)  <  u20) 

( 

p4l  =  P4(ullXl,tl+dl); 
p31  =P5(ulljU,tl,dl); 
p21  =P5(ullXl,tl-dl,dl); 
pll  =  P3(ullXl,tl-dl); 
p40  =  Pl(ulOXl.tl+dl); 
p30  =  P6(ul0Xl,tl,dl); 
p20  =  P6(ulOJLl,tl-dl,dl); 
plO  =  P7(ulOXl,tl-dl); 

q41  =Pl(u21X2,t2+d2); 
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q31  =  P2(u21,L2,t2,d2); 
q21  =  P5(u21,L2,t2-d2,d2); 
ql  1  =  P3(u21£2,t2-d2); 
q40  =  Pl(u20JL2,t2+d2); 
q30  =  P6(u20,L2,t2,d2); 
q20  =  P2(u20,L2,t2-d2,d2); 
qlO  =  P3(u20X-2,t2-d2); 

retum(  U(p4 1  ,p3 1  ,p2 1  ,p  1 1  ,p40,p30,p20,p  10,q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,q  10, 

k)); 


if  ((tl+dl)  <  ull  &&  (tl-dl)  <  ulO  &&  (t2+d2)  <  u21  &&  (t2-d2)  >=  u20) 

{ 

p41  =P4(ull,Ll,tl+dl); 
p31  =P5(ullJ-l,tltdl); 
p21=P5(ulU-l,tl-dl,dl); 
pi  1  =  P3(ul  1,L1, tl-dl); 
p40  =  Pl(ulOJLl.tl+dl); 
p30  =  P6(ulO,Ll,tl,dl); 
p20  =  P2(ul0,Ll,tl-dl,dl); 
plO  =  P3(ul0JLl,il-dl); 


q41  =  P4(u21J-2,t2+d2); 
q31  =  P5(u21  JL2,t2,d2); 
q21«P5(u21,L2,t2-d2,d2); 
qll  =  P3(u21,L2,t2-d2); 
q40  =  Pl(u20J-2,t2+d2); 
q30  =  P6(u20,L2.t2,d2); 
q20  =  P6(u20^2,t2-d2,d2); 
qlO  =  P7(u20,L2,t2-d2); 

rctum(  U(p4 1  ,p3 1  ,p2 1  ,p  1 1  ,p40,p30,p20,p  10,q4 1  ,q3 1  ,q2 1  ,q  1 1  tq40,q30,q20,q  10, 
k)); 

} 

if((tl+dl)<ull  &&  (tl-dl)<  ulO  &&  (t2+d2)<u21  &&  (t2-d2)<u20) 

{ 

p41  =  P4(ull,Ll, tl+dl); 
p31  =  P5(ull,Ll,tl,dl); 
p21  =  P5(ullJLl,tl-dl,dl); 
pll  =  P3(ull  J-l.tl-dl); 
p40  =  Pl(ul0,Ll,tl+dl); 
p30  =  P6(ul0,Ll.tl.dl); 
p20  =  P2(ul0,Ll,tl-dl,dl); 
plO  =  P3(ul0J-l,tl-dl); 
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q41 »  P4(u21,L2,t2+d2); 
q31  =  P5(u21,L2,t2,d2); 
q21  *  P5(u214-2,t2-d2,d2); 
ql  1  =  P3(u21J-2,t2-d2); 
q40  =  Pl(u20,L2,t2+d2); 
q30  *  P6(u20,L2,t2,d2); 
q20  =  P2(u20,L2,t2-d2,d2); 
qlO  =  P3(u20JL2,t2-d2); 

retum(  U(p4  l,p3  l,p21  ,pl  I,p40,p30,p20,pl0,q41  ,q3 1  ,q21  ,ql  1  ,q40,q30,q20,ql0, 

k)); 

} 

if((tl+dl)<ull  &&  (tl-dl)  <  ulO  &&  (t2+d2)>=u2l  &&  (t2-d2>=u20) 

{ 

p4I  =  P4(ullJ-l,tl+dl); 
p31  =  P5(ulUl,tl,dl); 
p21  =P5(ullJ-l, tl-dl  jdl); 
pll  =  P3(ul  1  JLl.tl-dl); 
p40  =  Pl(ulOJ-l,tl+dl); 
p30  =  P6(ulO,Ll,tl,dl); 
p20  =  P2(ulOJ-l,tl-dl,dl); 
pIO  =  P3(ulO,Ll,il-dl); 

q41  *  Pl(u21JL2,t2+d2); 
q31  =  P2(u21  JL2,t2,d2); 
q21  =  P5(u21JL2,t2-d2,d2); 
qll  =  P3(u21JL2,t2-d2); 
q40  =  Pl(u20JL2,t2+d2); 
q30  =  P6(u20^L2,t2,d2); 
q20  =  P6(u20,L2,t2-d2,d2); 
qlO  =  P7(u20JL2,t2-d2); 

retum(  U(p41,p3  l,p21,pl  1  ,p40,p30,p20,pl0,q4  l,q3 1  ,q21  ,ql  1  tq40,q30,q20,ql0, 

k)); 

} 

if((tl+dl)<ull  &&  (tl-dl)  <  ulO  &&  (t2+d2)>=u21  &&  (t2-d2)<u20  ) 

{ 

p41  =  P4(ullJ-l,tl+dl); 
p31  =  P5(ullJLl,tl,dl); 
p21  =  P5(ullJLl,tl-dl,dl); 
pll  =P3(ullJLl,tl-dl); 
p40  =  Pl(ui0,Ll,tl+dl); 
p30  =  P6(ulOJLl,tl,dl); 
p20  =  P2(ul0JLI,tl-dl,dl); 
plO  =  P3(ulOJLl,tl-dl); 

q41  =  Pl(u21JL2,t2+d2); 
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q31  =  P2(u21,L2,t2,d2); 
q21  =  P5(u21rL2,t2-d2,d2); 
ql  1  =  P3(u21,L2,t2-d2); 
q40  =  Pl(u204-2,t2+d2); 
q30  =  P6(u20J-2,t2,d2); 
q20  =  P2(u20J_2,t2-d2,d2); 
qlO  =  P3(u20,L2,t2-d2); 

retum(  U(p4 1  ,p3 1  ,p21  ,p  1 1  ,p40,p30,p20,pl0,q4 1  ,q3 1  ,q2 1  ,q  1 1  tq40,q30,q20,q  10, 
k)); 


N ****** **** **** / 


if  ((tl+dl)  >=  ul  1  &&  (tl-dl)  >=  ulO  &&  (t2+d2)  <  u21  &&  (t2-d2)  >=  u20) 

{ 

p4l  =  Pl(ull,Ll,tl+dl); 
p31  =  P2(ulU.l,tl,dl); 
p21  =  P5(ullJ-l,tl-dl,dl); 
pll  =P3(ullJLl,tl-dl); 
p40  =  Pl(ul0JLl, tl+dl); 
p30  =  P6(ulOJLl,tl  ,dl); 
p20  =  P6(ulOJLl,tI-dI,dI); 
plO  =  P7(ul0JLl,tl-dl); 


q41  =  P4(u21X2.t2+d2); 
q31  =  P5(u21,L2,t2,d2); 
q21  =  P5(u2 1  JL2,t2-d2,d2); 
qll  =P3(u21,L2,t2-d2); 
q40  =  P 1  (u20JL2,t2+d2); 
q30  -  P6(u20JL2,t2,d2); 
q20  =  P6(u20JL2,t2-d2,d2); 
qlO  =  P7(u20JL2,t2-d2); 

retum(  U(p4 1  ,p3 1  ,p21  ,p  1 1  ,p40,p30,p20,p  1 0,q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,q  10, 
k)); 

) 

if((tl+dl)>=ull  &&  (tl-dl)>=  ulO  &&  (£2+d2)<u21  &&  (t2-d2)<u20 ) 

{ 

p41  =  Pl(ullJLl,tl+dl); 
p31  =  P2(ul  lJLl.tl.dl); 
p21  =  P5(ullJLl,tl-dl,dl); 
pi  l  =  P3(ullJLl,  tl-dl); 
p40  =  Pl(ulOJLI,tl+dl); 
p30  =  P6(ul0,Ll,tl,dl); 
p20  =  P6(ul0,Ll, tl-dl, dl); 
pl0  =  P7(ul0Xl, tl-dl); 


93 


q41  =  P4(u21X2,l2+d2); 
q31  =  P5(u21X2,t2,d2); 
q21  =  P5(u21JL2,t2-d2,d2); 
qll  =  P3(u21JL2,t2-d2); 
q40  =  Pl(u20X2,t2+d2); 
q30  =  P6(ii2GjL2,t2,d2); 
q20  =  P2(u20JL2,t2-d2,d2); 
qlO  =  P3(u20X2,l2-d2); 

retum(  U(p4 1  ,p3 1  ,p21  ,p  1 1  ,p40,p30,p20,p  10, q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,ql0, 

k )); 

} 

if((tl+dl>=ull  &&  (tl-dl)  >=  ulO  &&  (t2+d2)>=u21  &&  (l2-d2)>=u20) 

{ 

p41  =  Pl(ul  1  JLl,tl+dl); 
p31  =  P2(ullJLl,tl,dl); 
p21  *  P5(ullJLl»tl-dl,dl); 
pll=P3(ullXl,tl-dl); 
p40  =  Pl(ul0Xl.tl+dl); 
p30  =  P6(ul0Xl,tl.dl); 
p20  =  P6(ulOXl,tl-dl,dl); 
plO  =  P7(ulOJLl,tl-dl); 

q41  a  Pl(u21JL2,t2+d2); 
q31  =  P2(u21X2,t2,d2); 
q21  =  P5(u21X2,t2-d2,d2); 
ql  I  =  P3(u21X2,t2-d2); 
q40  =  Pl(u20X2,t2+d2); 
q30  =  P6(u20X2,t2,d2); 
q20  =  P6(u20X2.t2-d2,d2); 
qlO  =  P7(u20X2,t2-d2); 

retum(  U(p4 1  ,p3 1  ,p21  ,p  1 1  ,p40,p30,p20,p  10, q4 1  ,q3 1  ,q2 1  ,ql  1  ,q40,q30,q20,q  10, 

k)); 

} 

if((tl+dl)  >=  ul  1  &&  (tl-dl)  >=  ulO  &&  (t2+d2)>=u21  &&  (t2-d2)<u20 ) 

{ 

p41  =  Pl(ullXl,tl+dl); 
p31  =  P2(ullXl,tl,dl); 
p21  =  P5(ul  1X1, tl-dl, dl); 
pll  =P3(ullXl,tl-dl); 
p40  =  Pl(ul0Xl,tl+dl); 
p30  =  P6(ulOXl,ll,dl); 
p20  =  P6(ul0Xl,tl-dl,dl); 
plO  =  P7(ulOXl,tl-dl); 

q41  =  Pl(u21X2,l2+d2); 
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q31  -  P2(u21X2,t2,d2); 
q21  *  P5(u21X2.t2-d2,d2); 
qll  =  P3(u21X2,t2-d2); 

•q40  =  Pl(u20X2.t2+d2); 
q30  =  P6(u20X2,t2,d2); 
q20  =  P2(u20X2,t2-d2,d2); 
qlO  =  P3(u20X2,t2-d2); 

retum(  U(p4  l,p3  l,p21,pl  1  ,p40,p30,p20,pl0,q41  ,q3 1  ,q21  ,ql  1  ,q40,q30,q20,ql0, 

k)); 

j y************* ****** *************** ************ *********************^ 

if  ((tl+dl)  >=  ul  1  &&  (tl-dl)  <  ulO  &&  (t2+d2)  <  u21  &&  (t2-d2)  >=  u20) 

{ 

p41=Pl(ullXl,tl+dl); 
p31  =  P2(ullJLl,tl,dl); 
p21»P5(ullJLl,tl-dl,dl); 
pll  =  P3(ullXl,tl-dl); 
p40  =  Pl(ulOJ^l,tl+dl); 
p30  =  P6(ulOJLl,tl,dl); 
p20  =  P2(ul0J.l,tl-dl41); 
plO  =  P3(ulOJ-l,tl-dl); 

q41=P4(u21X2,t2+d2); 
q31«P5(u21JL2,t2,d2); 
q21  =  P5(u21X2,t2-d2,d2); 
ql  1  =  P3(u21  JL2,t2-d2); 
q40  =  Pl(u20X2,t2+d2); 
q30  -  P6(u20JL2,t2,d2); 
q20  =  P6(u20JL2,t2-d2,d2); 
qlO  =  P7(u20X2,t2-d2); 

retum(  U(p4 1  ,p3 1  ,p21  ,p  1 1  ,p40,p30,p20,p  10, q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,q  10, 

k)); 

} 

if((tl+dl)>=ul  1  &&  (tl-dl)  <  ulO  &&  (t2+d2)<u21  &&  (t2-d2>cu20 ) 

( 

p41  =  Pl(ul  1X1, tl+dl); 
p31  =  P2(ullXl,tl,dl); 
p21  =  P5(ullXl,tl-dl,dl); 
pi  l  =  P3(ul  1X1,  tl-dl); 
p40  =  Pl(ul0XUl+dl); 
p30  =  P6(ul0Xl,tl.dl); 
p20  =  P2(ulOXl,tl-dl,dl); 
plO  =  P3(ulOXl,tl-dl); 
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q41=P4(u21,L2.t2+d2); 
q31  =  P5(u2U_2.t2,d2); 
q21  =  P5(u21X2,t2-d2,d2); 
qll  =  P3(u21,L2,t2-d2); 
q40  =  Pl(u20J-2,t2+d2); 
q30  =  P6(u20,L2,t2,d2); 
q20  =  P2(u20J-2,t2-d2,d2); 
qlO  =  P3(u20J-2,t2-d2); 

retum(  U(p4 1  ,p3  l,p21  ,p  1 1  ,p40,p30,p20,pl0,q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,ql0, 

k»; 

} 

if((tl+dl)>=ull  &&  (tl-dl)  <  ulO  &&  (t2+d2)>=u21  &&  (t2-d2)>=u20) 

{ 

p41=Pl(ullJ-l,tl+dl); 
p31  =  P2(ull,Ll,tl,dl); 
p2I  =  P5(ullJLl,tl-dl,dl); 
pH  =  P3(ul  IX- 1,  tl-dl); 
p40  =  Pl(ulO,Ll,tl+dl); 
p30  =  P6(ul0Xl,tl,dl); 
p20  =  P2(ul0Xl,tl-dl,dl); 
plO  =  P3(ulO,Ll,tl-dl); 

q41  =  Pl(u21JL2,t2+d2); 
q31  =  P2(u21X2,t2,d2); 
q21  =  P5(u21,L2,t2-d2,d2); 
qll  =  P3(u21,L2,t2-d2); 
q40  =  Pl(u20X2.t2+d2); 
q30  =  P6(u20JL2,l2,d2); 
q20  =  P6(u20,L2,t2-d2,d2); 
qlO  =  P7(u20,L2,t2-d2); 

retum(  U(p4 1  ,p3 1  ,p21  ,p  1 1  ,p40,p30,p20,p  10,q4 1  ,q3 1  ,q2 1  ,q  1 1  ,q40,q30,q20,q  10, 

k)); 

} 

if((tl+dl)  >=  ul  1  &&  (tl-dl)  <  ulO  &.&  (t2+d2)>=u21  &&  (t2-d2)<u20 ) 

{ 

p4I  =  Pl(ullrLl,tI+dl); 
p31  =  P2(u  1 1  1  ,t  1  ,d  1 ); 
p21  =  P5(ullXl.tl-dl.dl); 
pll  =P3(ullJLl,ll-dl); 
p40  =  Pl(ulO,Ll,tl+dl); 
p30  =  P6(ulO,Ll,tl,dl); 
p20  =  P2(ulO,Ll,tl-dl,dl); 
plO  =  P3(ulO,Ll, tl-dl); 

q41  =  Pl(u21,L2,t2+d2), 
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q31  =  P2(u21,L2,t2,d2); 
q21  =  P5(u21^2,t2-d2,d2); 
qll  =  P3(u21,L2,t2-d2); 
q40  =  Pl(u20X.2,t2+d2); 
q30  =  P6(u20JL2,t2,d2); 
q20  =  P2(u20J_2,t2-d2,d2); 
qlO  »  P3(u20X2,t2-d2); 

rctum(  U(p4 1  ,p3 1  ,p21  ,p  1 1  ,p40,p30,p20,p  1 0,q4 1  ,q3 1  ,q2 1  ,q  1 1  tq40,q30,q20,q  10, 

k)); 

}  ) 


/**********•**********•********************************/ 


/•calculate  constants  in  the  G  function  •/ double  Cons  10  { 
double  a,b; 
double  A3; 
double  i; 

A  =  (beta)  /  (1  -  alpha); 

B  =  (1  *  beta)  /  alpha; 
a  =  log(A); 
b  =  log(B); 

i  =  (a*A*(B-l)  +  b*B*(l-  A))  /  (B  -  A); 
retum(i); ) 

double  ConsOO  { 
double  a,b; 
double  A3; 
double  i; 

A  =  (beta)  /  (1  -  alpha); 

B  =  (1  -  beta) /alpha; 
a  =  log(A); 
b  =  log(B); 

i  =  (a  *  (B  -  1)  +  b  *  (1  -  A))  /  (B  -  A); 
retum(i);  ) 
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APPENDIX  C 


Program  for  Evaluation  of  G  Function  of  Two-level  Quantizer  System 


This  program  is  to  generate  double  exponential  random  environment,  and  it  is 
then  be  used  to  simulate  the  sequential  detection  system  with  two  and  four  level 
quantization.  The  number  of  data  that  have  to  be  collected  to  determine  which  signal 
has  been  received  is  calculated.  This  program  makes  comperison  between  2  and  4 
level  quatizers. 

Since  these  two  sensors  are  independent,  the  random  generator  is  started  twice 
in  order  to  make  the  environment  of  sensor  one  independent  from  the  one  of  sensor 
two. 
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#define  RND(seed)  rand(seed)  /  32767.  #define  RND(weed)  rand(weed)  /  32767. 
#define  max  500  #include  <math.h>  #include  <stdio.h>  int  randO;  void  srand(); 
unsigned  int  seed;  unsigned  int  weed; 

main()  { 
float  yl,y2; 
float  lamda.lamdal; 
float  u0,ul,U0,Ul; 

float  al  I,al2,a01,a02;  /*  quantized  areas  for  2  level  quantizer  */ 

float  arll,arl2,arl3,arl4,ar01,ai02,ar03,ar04;  /*  quantized  areas  for  4  level*/ 

float  All, A12,A01,A02; 

float  Arl  l,Arl2,Arl3,Arl4,Ai01,Ar02,Ai03,Ar04; 
float  tl,t2; 

float  Tl,T2,Deltal,  Delta2; 

float  fl,  fO; 

float  Fl,  FO; 

float  ffl,  ffO; 

float  FF1,  FFO; 

float  A,B; 

float  missl,miss2; 

float  testl,test2; 

float  airivlO,  arriv2(); 

float  Areal 0,  Area2(),Area3(),  Area4(),  Area5(),  Area6(),  Area70; 

float  alpha,  beta; 

int  ij,k,m,n; 

int  okl,ok2; 

int  Nl[max+1]; 

int  N2[mar+lj; 

FILE  *fpl,  *fp2,  *fp3,  *fp4; 

/*  initial  the  random  number  for  sensor  1  */ 

seed  =  time(&seed); 

srand(seed);  /*Given  the  constants  of  double  exponential  function*/ 
lamda  =  0.25; 
lamdal  =  0.125; 
uO  =  -3.0; 
ul  =  3.0; 

U0  =  -5.0; 

Ul  =  5.0;  /‘define  the  optimal  quantizers  of  for  sensor  1  and  2  */ 
tl  =  0.0; 
t2  =  0.0; 

Tl  =  0.0; 

T2  =  0.0; 

Delta  1  =  2.0; 

Delta2  =  3.3; 
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alpha  =  0.05; 

beta  =  0.01;  /*calculate  the  threshold  of  sequential  test*/ 

A  =  beta  /  ( 1  -  alpha ); 

B  =  ( 1  -  beta /  alpha ; 
fpl  =  fopen("numlfig4","w"); 
fp2  =  fopen("num2fig4","w"); 
fp3  =  fopen("mislfig4","w"); 

fp4  =  fopen("mis2fig4","w");  /*find  the  probabilities  of  each  quantized  level  for  a 
given  */  /*  quantizer  */ 

all  =  Area2(lamda,ul,tl); 
al2  =  Area3(lamda,ul,tl); 
aOl  =  Area4(lamda,u0,tl); 
a02  =  Areal(lamda,uO,tl); 

All  =  Area2(lamdal,Ul,t2); 

A12  =  Area3(lamdal,Ul,t2); 

AOl  =  Area4(lamdal,U0,t2); 

A02  =  Areal(lamdal,U0,t2); 
if  (( T1  +  Deltal)  >=ul) 

{ 

aril  =  Area2(lamda,ul,Tl-Deltal); 
arl2  =  Area50amda,ul,Tl-Deltal,Tl); 
arl3  =  Area7(lamda,ul,Tl,Deltal+Tl); 
arl4  =  Areal(lamda,ul,Tl+Deltal); 
arOl  =  Area20amda,uO,Tl -Deltal); 
ar02  =  Area7(lamda,uO,Tl-Deltal,Tl); 
ar03  =  Area6(lamda,uO,Tl,Tl+Deltal); 
ai04  =  Area  1  (lamda,uO,T  1  +Del  ta  1 ) ; 

} 

if  (( T1  +  Deltal  )<ul ) 

{ 

aril  =  Area2(lamda,ul,Tl -Deltal); 
arl2  =  Area5(lamda,ul,Tl-Deltal,Tl); 
arl3  =  Area5(lamda,ul,Tl,Deltal+Tl); 
arl4  =  Area3(lamda,u  1  ,T  1  +Delta  1 ); 
arOl  =  Area4(lamda,uO,Tl-Deltal); 
ar02  =  Area6(lamda,uO,Tl-Deltal,Tl); 
ar03  =  Area6(lamda,uO,T  1  ,T  1  +Delta  1 ); 
ar04  =  Areal(lamda,uO,Tl+Deltal); 

} 

if(  (T2  +  Delta2)  >=  U1 ) 

{ 

Aril  =  Area2(lamdal,Ul,T2-Delta2); 

Arl2  =  Area5  (lamda  1  ,U  1  ,T2-Del  ta2,T2); 

Arl3  =  Area7  (lamda  1  ,U  1  ,T2,T2+Delta2); 

Arl4  =  Areal(lamdal,Ul,T2+Delta2); 
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ArOl  =  Area2(lamdal  ,U0,T2-Delta2); 

Ar02  =s  Arca7(lamdal,UO,T2-Dclta2,T2); 

Ar03  =  Area6(Iamdal,U0,T2,T2+DeIta2); 

Ar04  =  Area  1  (lamda  1  ,U0,T2+Delta2); 

} 

if(  (T2  +  Delta2)  <  U1 ) 

{ 

Aril  =  Area2(lamdal,Ul,T2-Delta2); 

Arl2  =  Area5(lamda  1  ,U  1  ,T2-Delta2,T2); 

Arl3  =  Area5(lamda  1  ,U  1  ,T2,T2+Delta2); 

Arl4  =  Area3(lamda  1  ,U  1  ,T2+Delta2); 

ArOl  =  Area4(lamda  1  ,U0,T2-Delta2); 

Ar02  a  Area6(lamdal  ,U0,T2-Delta2,T2); 

Ar03  =  Anea6(lamdal  ,U0,T2,T2+Dclta2); 

Ai04  a  Area  1  (lamda  1  ,UO,T2+Delta2); 

)  l*  initial  the  random  number  for  sensor  2  */ 

weed  =  time(&weed); 
srand(weed); 

m  =  0; 
n  =  0; 

for(i=l;  i<a  max;  i++) 

{ 

j  —  0; 
k  =  0; 
okl  =0; 
ok2  =  0; 
lestl  =  1.0; 
test2  a  1.0; 

while((okl  =  0)  II  (ok2  =  0)) 

{ 

yl  =  arrivl(lamda,ul); 
y2  =  arriv2(lamdal,Ul); 

f*  calculate  the  number  of  collected  data  from  two-level  quantizer  */ 
if(okl  =0) 

{ 


if(yl  >=  tl) 

{ 

fl  =  al2; 
f0  =  a02; 

} 

else  if  (yl  <=  tl) 


1DJ 


fl  =  all; 
fO  =  aOl; 

} 

if(y2  >=  (2) 

[ 

Fl  =  A 12; 

FO  =  A02; 

} 

else  if  (y2<=  1 2) 

{ 

Fl  -All; 

FO  =  AOl; 

} 

testl  =  ((fl  *  Fl)  /  (fO  *  FO))*testl; 
if  ((test  1  >=  B)  II  (testl  <=  A)) 
okl  =  1; 
if(testl  <=  A) 

{ 

m++; 

} 

j++; 

}  /*  calculate  the  number  of  collected  data  from  four-level  quantizer  */ 
if(ok2  =  0) 

{ 

if(yl  <  T1  -  Deltal) 

{ 

ffl  =  aril; 
ffO  =  aiOl; 

} 

if(y  1  >=  (T1  -  Deltal)  &&  yl  <  Tl) 

{ 

ffl  =  arl2; 
ffO  =  ar02; 

} 

if(yl  >=  Tl  &&  yl  <  (Tl  +  Deltal)) 

{ 

ffl  =  arl3; 
ffO  =  ai03; 

} 

if(yl  >=  (Tl  +  Deltal)) 

{ 

ffl  =  arl4; 
ffO  =  ar04; 

} 

if(y2  <  (T2  -  Delta2)) 
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{ 

FF1  =  Aril; 

FFO  =  ArOl; 

} 

if(y2  >=  (T2  -  Delta2)  &&  y2  <  T2) 

{ 

FFl  =  Arl2; 

FFO  =  Ai02; 

} 

if(y2  >=  T2  &&  y2  <  (T2  +  Delta2)) 

{ 

FFl  =  Arl3; 

FFO  =  Ar03: 

} 

if(y2  >=  (T2  +  Delta2)) 

{ 

FFl  =  Arl4; 

FFO  =  Ar04; 

} 

test2  =  ((ffl  *  FFl)  /  (ffO  *  FF0))*test2; 
if((test2  >=  B)  II  (test2  <=  A)) 
ok2  =  1; 
if(test2  <=  A) 
n++; 
k++; 

} 

) 

Nl[i]  =  j; 

N2[i]  =  k; 

fprintf(fpl,"%d  %dO,i,Nl[i]); 
fprintf(fp2,"%d  %dO,i,N2[i]); 

} 

fclose(fpl); 

fclose(fp2); 

fclose(fp3); 

fclose(fp4);  }  /*  generate  the  enviomment  HI  for  sensor  1 
float  1,  u; 


float  z,r; 
int  k,y; 
double  x; 

x  =  RND(seed); 
if  (x  <  0.5  &&  x  >  0.0) 
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*/  float  arrivl(l,  u  ) 


{ 

y  =  i; 

z  =  u  +  (y  *  log(2  *  x))  / 1 ; 
retum(z); 

) 

else  if  ( x  >=  0.5  &&  x  <  1.0) 

{ 

y  =  -i; 

z  =  u  +  (y  *  log  (2*(l-x)))/l; 
retum(z); 

}  }  /*  generate  the  enviomment  HI  for  sensor  2  */  float  arriv2(l,  u  )  float  1,  u; 


{ 

float  z,r, 
int  k,y; 
double  x; 

x  =  RND(weed); 
if  (x  <  0.5  &&  x  >  0.0) 

{ 

y  =  l; 

z  =  u  +  (y  *  log(2  *  x))  / 1 ; 
retum(z); 

} 

else  if  (x  >=  0.5  &&  x  <  1.0 ) 

{ 

y  =  -i; 

Z  =  U  +  y  *  (log  ( 2  *  (  1  -  x )))  /I ; 
retum(z); 

)} 

I*  calculate  the  probabilies  of  different  quantized  region  for  */  /*  two  and  four-level 
quantizer  */ 

float  Areal  (l,u,x)  float  l,u,x;  ( 

retum(0.5  *  exp(-l  *  (x  -  u)));  }  float  Area2(l,u,x)  float  I,u,x;  { 
retum(0.5  *  exp(-l  *  (u  -  x)));  }  float  Area3(l,u,x)  float  l,u,x;  { 
retum(l  -  0.5  *  exp(-l  *  (u  -  x)));  }  float  Area4(l,upc)  float  I,upc;  { 
retum(l  -  0.5  *  exp(-l  *  (x  -  u)));  )  float  Area5(l,u,x,y)  float  l,u,x,y;  { 
retum(0.5  *  exp(-l  *  (u  -  y))  -  0.5  *  exp(-l  *  (u  -  x)));  }  float  Area6(l,u,x,y)  float 
l,u,x,y;  { 

retum(0.5  *  exp(-l  *  (x  -  u))  -  0.5  *  exp(-l  *  (y  -  u)));  )  float  Area7(l,u,x,y)  float 
l,u,x,y;  { 

retum(l  -  0.5  *  exp(-l  *  (u  -  x))  -  0.5  *  exp(-l  *  (y  -  u)));  } 
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MISSION 

OF 

ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search t,  development,  test,  and  technology  transition  in  support  of  Air 

O 

Force  Command,  Control,  Communications  and  Intelligence  (C  l)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
ESD  elements  to  perform  effective  acquisition  of  C3I  systems.  In  addition, 
Rome  Laboratory’s  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/ maintainability  and  testability. 


